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1  FISCAL  YEAR  201 5  TECHNICAL  ACHIEVEMENTS 


1.1  Research  Objectives 

Research  objectives  for  this  fiscal  year  include: 

•  Quantify  and  develop  scaling  rules  for  mass  transport  and  heat  release  distribution  in  ultra¬ 
compact  combustors 

•  Characterization  of  physical  and  chemical  processes  controlling  stability  of  high-g 
combustors 

•  Assess  limitations  of  scaling  high-g  combustors  to  larger  sizes 

•  Quantify  distribution  and  interaction  of  fuel  droplets  in  high  centrifugal  force  environments 

•  Flame  distribution  and  heat  flux  interrogation  during  reacting  film  cooling 

•  Continue  experimental  and  computational  investigations  of  fundamental  combustion 
phenomena. 

•  Explore  the  details  of  turbulent  mixing  and  turbulence-chemistry  interactions  through  studies 
of  fundamental  turbulent  events. 

•  Investigate  the  physics  and  chemistry  of  combustion  processes  responsible  for  the  formation 
of  particulates. 

•  Characterize  combustion  ignition  phenomena  and  development  of  novel  ignition  schemes. 

•  Develop,  demonstrate,  and  apply  advanced  laser-based  and  other  optical  diagnostic 
techniques  for  the  experimental  characterization  of  key  combustion  phenomena. 

•  Develop,  demonstrate,  and  apply  advanced  modeling  and  simulation  (M&S)  techniques  for 
the  computational  characterization  of  key  combustion  phenomena. 

•  Continue  efforts  to  integrate  increasingly  sophisticated  chemistry  and  kinetics  models  with 
fluid  dynamics  to  achieve  superior  combustion  codes. 

•  Achieve  a  validated  predictive  model  of  turbulent  mixing  and  turbulence-chemistry 
interactions. 

1.2  Background  of  Research 

During  the  past  fiscal  year,  researchers  with  the  Air  Force  Research  Laboratory’s  Combustion 
Branch  made  substantial  progress  in  numerous  areas  including:  1)  ultra-compact  combustors,  2) 
inter-turbine  burner  concepts,  3)  bluff-body  stabilized  turbulent  flames,  4)  well-stirred  reactors 
for  chemical  kinetics,  and  5)  detonation-stabilized  turbulent  flames. 

The  approaches  used  to  meet  the  objectives  of  this  program  are  as  follows.  First,  experiments  are 
performed  in  state-of-the-art  atmospheric  and  high  pressure/temperature  facilities.  A  suite  of 
diagnostics  including  pressure  and  temperature  measurements,  gas  sampling  and  analysis,  visual 
imaging,  and  optical  (laser)  measurements,  are  used  to  understand  flame  phenomenon  and 
characterize  the  performance  of  combustion  systems.  Fundamental  flame  studies  provide  insights 
into  fluid  mechanic  and  chemistry  effects  within  combustors.  Larger  scale  experiments  are  used 
to  validate  combustor  components  and  to  bridge  results  from  the  fundamental  flame  studies  to 
engine  conditions.  The  second  approach  for  meeting  the  research  objectives  is  to  perform 
modeling  and  simulations  in  conjunction  with  the  experiments.  The  calculations  provide  further 
insights  into  combustion  phenomenon  and  explore  concepts  which  are  not  viable  experimentally. 
Computational  work  includes  two-dimensional  calculations  with  full  combustion  chemistry  and 
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three-dimensional  simulations  which  capture  fluid  mechanic  effects.  Parametric  studies  are  used 
to  explore  new  design  space  and  improve  experiments.  The  third  aspect  of  meeting  the  research 
objectives  is  using  SBIR  programs  to  support  developing  new  technologies  and  transitioning  this 
work  to  the  warfighter.  In-house  testing  of  SBIR  products  supports  validation  of  innovative 
concepts  and  leverages  the  high  and  low  pressure  testing  capabilities  within  the  turbine  engine 
division. 

The  "Highly  Loaded,  Low  Emissions  Combustion  Systems"  program  encompasses  a  range  of 
projects  from  more  fundamental  to  applied  tests.  The  combined  approach  of  the  UCC  efforts  and 
span  of  research  allows  a  pathway  for  aerodynamic  integration  of  compression,  combustion,  and 
turbine  airflows.  This  combination  integrates  the  compressor  exit  guide  vane  along  with  the  first 
stage  turbine  nozzle  with  a  more  compact  combustor.  The  combination  reduces  airflow  losses  by 
decreasing  flow  angular  momentum  changes,  resulting  in  a  reduced  pressure  drop  and  greater 
engine  performance  in  a  decreased  volume.  These  compact  combustors  also  enable  their  use 
between  the  high  pressure  and  low  pressure  turbines,  resulting  in  a  revolutionary,  reheat  cycle 
propulsion  system  called  the  ITB.  The  UCC  has  the  potential  for  a  2.4  percent  decrease  in  engine 
weight,  0.8%  reduction  in  thrust  specific  fuel  consumption,  and,  with  an  ITB,  a  4,500  hp 
potential  for  additional  power  extraction  at  cruise  in  comparison  to  a  Gen  V  fighter  engine.  This 
approach  can  revolutionize  engine  design  concepts  for  Gen  VI  fighter  engines. 

The  AFRL  team  was  the  first  in  the  world  to  accomplish  several  groundbreaking  milestones 
towards  achieving  a  practical  UCC.  Multiple  tests  of  a  UCC  were  performed  at  engine  relevant 
conditions.  The  technology  demonstrated  that  a  single  cavity  can  be  used  for  the  combustor 
while  maintaining  appropriate  heat  release  and  temperature  profiles  for  engines.  This 
advancement  will  reduce  the  number  of  components  needed  for  a  combustor,  thus  decreasing  the 
cost  and  complexity  of  the  system.  Evaluation  of  key  fluid  mechanic  interactions  between  the 
turbine  components  and  the  combustor  were  performed  and  used  to  improve  design  rules  for  the 
system  level  approach. 

The  team  developed  a  series  of  novel  design  and  analysis  tools  to  enable  development  and 
facilitate  transitioning  technology  to  OEMs.  They  designed  and  implemented  a  data  processing 
technique  to  quantify  visible  observations  from  high-speed  videos  and  provide  correlations  with 
measured  combustor  performance.  Further,  the  team  acquired  planar  laser-induced  fluorescence 
(PLIF)  measurements  in  the  combustor;  some  of  only  a  few  sets  of  PLIF  measurements  collected 
at  engine  relevant  conditions.  The  PLIF  technique  spatially  measures  flame  location  and  heat 
release  in  the  combustor,  enabling  validation  of  modeling  codes  using  data  collected  non- 
intrusively  in  the  flame  zone.  Finally,  they  analyzed  a  novel  approach  to  reducing  the  number  of 
components  in  the  compressor  and  improving  the  UCC  system  for  the  first  time.  Scientific 
insights  gained  regarding  cooling  turbine  components  within  a  combustor,  fuel  droplet 
distribution,  and  interaction  of  flames  in  unconventional  system  will  reduce  the  cost  associated 
with  developing  and  transitioning  the  UCC  concept. 
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2 


RESEARCH  REVIEW 


2.1  Studies  on  the  Effects  of  Centrifugal  Force  on  the  Performance  of  a  High-G 
Ultra  Compact  Combustor 

2.1.1  Introduction 

There  is  a  need  to  develop  advanced  combustors  for  aircraft  gas  turbine  engines  with  improved 
thrust-to-weight  ratio,  reduced  gaseous  emissions,  and  lower  pressure  losses  and  specific  fuel 
consumption.  Such  combustors  must  preserve  low  maintenance  cost  and  high  durability.  In  order 
to  accomplish  this  goal  it  is  necessary  to  reduce  the  size  of  the  combustors.  Every  inch  of  a 
combustor  length  reduction  translates  to  significant  weight  savings  (this  is  highly  dependent  on 
the  engine  size  and  type)  [1].  Reductions  in  engine  length  and  weight  favorably  impact  engine 
thrust  to  weight  ratio  and  fuel  bum  performance.  A  combustor  with  shorter  length  and  smaller 
volume  equates  to  a  shorter  residence  time  compared  to  conventional  combustors.  Since  gaseous 
emissions  such  as  NOx  tend  to  form  in  long  residence  time  regions  at  elevated  temperature,  this 
is  an  advantage  because  these  emissions  would  have  less  time  to  form  within  the  combustor. 
However,  a  shorter  combustor  imposes  some  difficult  challenges  to  overcome.  A  shorter,  more 
compact  combustor  means  that  there  is  less  time  for  the  fuel  to  mix  with  combustor  air  and 
completely  bum  within  the  combustor.  Therefore,  faster  mixing  of  fuel  with  air  and  reactants 
with  products  as  well  as  increased  turbulent  flame  speeds  must  occur.  This  phenomenon  could 
potentially  lead  to  an  increase  in  pressure  losses  in  order  to  enhance  mixing. 

The  UCC  concept  combines  the  primary  and  the  dilution  zones  of  a  conventional  combustor  with 
high  pressure  turbine  inlet  guide  vanes  (IGVs)  to  shorten  the  length  of  the  combustor  compared 
to  a  conventional  design  with  separate  combustor  and  IGV  components.  There  are  two  general 
concepts  of  the  UCC,  viz.,  high-g  cavity  (HGC)  and  trapped- vortex  cavity  (TVC).  Both  concepts 
incorporate  a  circumferential  cavity  that  is  recessed  from  the  mainstream  flow.  Liquid  fuel  spray 
and  hot  air  are  injected  in  this  circumferential  cavity  (CC)  where  fuel  evaporates  and  bums  in  a 
fuel-rich  primary  zone.  Low  pressure  wakes  transport  the  products,  unbumed  fuel  and 
intermediate  species  from  the  cavity  to  the  mainstream  flow  to  complete  burning  in  a  secondary 
zone  at  fuel-lean  equivalence  ratios.  The  reacted  mixture  then  flows  through  the  IGVs,  which  act 
as  a  nozzle  to  increase  the  Mach  number  and  swirl  angle  for  entry  into  the  high  pressure  turbine 
rotor  stage  for  power  extraction.  The  HGC  and  TVC  differ  primarily  in  the  manner  in  which  fuel 
and  air  are  injected  in  the  recessed  cavity  [1,2].  For  the  UCC-TVC  concept,  air  is  axially  injected 
into  the  cavity  in  non-aligned  opposite  directions,  leading  to  a  cavity-stabilized  flame  vortex  with 
its  axis  tangentially  aligned  to  the  circumference  of  the  combustor.  Therefore,  there  is  no  actively 
forced  bulk  circumferential  swirl  flow  in  this  concept.  On  the  other  hand,  the  UCC-HGC  cavity 
air  is  injected  inwards  along  the  outer  wall  of  the  recessed  cavity  at  an  angle  to  the  tangent  of  the 
cavity  outer  wall  surface,  leading  to  a  bulk  circumferential  flow  within  the  cavity  around  the 
combustor  annulus.  This  bulk  swirl  in  the  cavity  generates  a  high  centrifugal  force  that  acts  as  a 
buoyant-like  force.  It  impulses  the  hot,  lower-density  reacted  gases  in  the  recessed  cavity  radially 
inwards,  while  pushing  the  colder,  higher-density  non-reacted  gases  radially  outward.  This 
allows  for  a  shorter  flame  that  oxidizes  fuel  at  a  “bubble”  flame  speed  [1-3]  faster  than  its 
corresponding  fuel-air  turbulent  flame  speed. 
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The  concept  of  the  HGC  was  conceived  from  the  experiments  of  Lewis  in  the  1970’s.  Lewis  and 
coworkers  [4-6]  used  a  cylindrical  tube  that  was  filled  with  a  premixed  mixture  of  fuel  and  air. 
Propane  and  hydrogen  were  used  in  these  experiments.  Like  a  propeller,  the  tube  was  rotated 
about  its  center,  and  once  it  had  reached  its  intended  rotational  speed,  an  igniter  positioned  at  one 
end  of  the  tube  sparked  a  flame  in  the  fuel-air  mixture.  The  centrifugal  force  induced  Rayleigh- 
Taylor  instability  by  negatively  stratifying  the  mixture  density  in  the  direction  of  this  force. 
Therefore,  the  centrifugal  force  is  measured  in  factors  of  the  Earth’s  gravitational  acceleration. 
Depending  on  the  calculated  g-load,  the  flame  propagated  towards  the  other  end  of  the  tube  at  a 
flame  propagation  velocity  equal  or  greater  than  the  mixture’s  corresponding  turbulent  flame 
speed.  They  claimed  that  at  low  centrifugal  forces  (less  than  500  g’s)  there  was  no  effect  on  the 
propane-air  flame  propagation  velocity.  However,  when  the  centrifugal  force  was  increased 
substantially  above  500  g,  they  claimed  that  the  turbulent  flame  speed  increased  by  up  to  a  factor 
of  four.  They  reasoned  that  the  flame  propagated  at  a  “bubble”  velocity  that  exceeded  its 
corresponding  turbulent  flame  speed.  The  sensitivity  of  the  flame  speed  of  hydrogen  to  g-loads 
was  impacted  by  the  equivalence  ratio  of  the  mixture.  The  flame  speed  of  stoichiometric 
hydrogen-air  was  independent  of  g-load  whereas  that  of  leaner  premixed  hydrogen-air  showed 
sensitivity  to  g-load  similar  to  stoichiometric  propane-air  premixed  flames,  ft  was  stated  that  at 
very  high  centrifugal  forces  of  nearly  3500g  the  “bubble”  velocity  peaked  and  then  began  to 
decrease  abruptly  until  the  flame  extinguished  due  to  very  high  stretch  rates. 

Nonetheless,  recent  numerical  simulations  have  not  corroborated  the  effect  of  centrifugal  force 
on  turbulence  flame  speed,  ft  is  vital  to  point  out  that  for  a  dynamic  system  the  flame 
propagation  velocity  is  not  equal  to  the  turbulent  flame  speed.  Flame  propagation  velocity  is  the 
observed  flame  velocity  from  a  reference  frame  and,  hence,  it  is  the  summation  of  the  local  flow 
velocity  and  the  turbulent  flame  speed.  Katta  and  coworkers  [1]  have  previously  investigated  this 
problem.  They  simulated  Lewis’  experiments  by  imposing  a  body  force  in  the  momentum 
equation  parallel  to  the  flow  mainstream  motion.  However,  their  simulations  also  differ  from 
those  from  Lewis  and  coworkers  [4-6]  because  the  tube  was  opened  at  one  end.  This  did  not 
allow  for  pressure  rise  that  actually  occurred  in  the  experiments.  However,  the  flame  propagation 
characteristics  compared  well  with  Lewis  experiments.  Recently,  Briones  et  al.  [2]  numerically 
examined  the  effect  of  centrifugal  force  on  flame  propagation  velocity  in  a  rotating  tube  that 
more  closely  mimics  Lewis  and  coworkers  experiments  [4-6],  Even  though  there  was  relatively 
good  comparison  between  measured  [4-6]  and  predicted  flame  propagation  velocity  as  function 
of  centrifugal  force,  Rayleigh-Taylor  instability  with  thermal  expansion  and  turbulent  flame 
speed  constitute  (in  descending  order)  the  flame  propagation  mechanism.  In  addition,  Katta  and 
coworkers  [1]  conducted  a  two-dimensional  simulation  where  they  imposed  high-g  forces  in  a 
JP-8  fuel-rich  recessed  cavity  for  a  UCC  application.  They  found  that  the  high-g  forces  had 
nearly  negligible  effect  on  their  two-dimensional  numerical  results. 

Therefore,  it  appears  that  the  main  reason  for  developing  a  UCC-HGC  with  expected  shorter 
flame  that  allows  for  a  smaller  combustor  such  as  that  previously  investigated  [1-3]  is  not 
fundamentally  founded  (i.e.,  the  centrifugal  force  may  have  no  effect  on  the  turbulent  flame 
speed).  Moreover,  the  UCC-HGC  was  also  conceptualized  on  the  implicit  premise  that  faster 
turbulent  flame  speed  would  translate  to  shorter  flames.  This  premise  also  appears  to  be 
erroneous  because  turbulent  flames  are  composed  of  multiple  premixed  and  nonpremixed 
combustion  regimes.  The  former  is  sensitive  to  turbulent  flame  speed  whereas  the  latter  is  not. 
Even  though  turbulent  flame  speeds  are  important  for  flame  stabilization  and  liftoff  there  is  no 
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indication  that  turbulent  flame  speeds  are  relevant  for  flame  length.  For  instance,  Zelina  and 
coworkers  [9-1 1]  have  claimed  that  the  UCC-HGC  can  achieve  flame  length  reductions  in  the 
order  of  50-70%  when  compared  to  conventional  combustors.  Nevertheless,  it  is  not  possible 
from  Refs.  [9-11]  to  discern  how  that  conclusion  was  reached  since  turbulent  flame  length 
measurements  are  not  discussed.  Despite  these  controversies  the  UCC-HGC  is  still  of  interest 
because  Zelina  and  coworkers  [9-11]  have  demonstrated  that  the  UCC-HGC  could  attain  higher 
combustion  efficiencies  for  a  wider  range  of  loading  parameter,  higher  heat  release  rate  and 
larger  static  stability  than  conventional  combustors,  while  maintaining  reasonable  NOx 
emissions. 

In  order  to  further  reach  a  consensus  on  the  effect  of  the  centrifugal  force  on  the  performance  of 
the  UCC-HGC,  we  numerically  investigate  the  effect  of  centrifugal  forces  on  the  turbulent  flame 
speed  and  flame  length  using  the  UCC-HGC.  For  this,  however,  the  previous  UCC-HGC  rig  [9- 
1 1]  is  modified  in  order  to  operate  at  more  realistic  conditions  than  it  did  before.  Zelina  and 
coworkers  operated  the  UCC-HGC  rig  [9-1 1]  at  fairly  low  operating  pressures  (~  4  atm),  low 
inlet  temperatures  (~550  K),  low  Mach  numbers  (~  0.1 -0.2)  and  high  centrifugal  forces  (~  3000 
g’s)  mainly  because  the  rig  was  not  cooled  and  was  smaller  than  a  practical  combustor.  In  the 
current  investigation  we  have  scaled  and  air-cooled  the  UCC-HGC  rig  so  that  it  actually 
corresponds  to  a  combustor  of  a  relevant  gas  turbine  engine  cycle.  Now,  the  UCC-HGC  rig  can 
operate  at  higher  pressures  (~19  atm),  higher  inlet  temperature  (~  800  K)  and  higher  Mach 
number  (~  0.5).  The  drawback  is  that  this  more  practical  UCC-HGC  would  operate  at  lower 
nominal  centrifugal  forces  up  to  approximately  1700  g’s.  In  the  current  investigation, 
nonetheless,  the  centrifugal  force  will  be  varied  by  changing  the  circumferential  cavity  (CC)  air 
driver  jets  angle  and  pressure  drop  (i.e.,  mass  flow  rate  or  cavity  loading),  across  the  driver 
holes. 

The  purpose  of  this  investigation  is  to  provide  further  insights  on  whether  or  not  centrifugal 
forces  enhance  the  combustion  performance  of  the  UCC-HGC.  Therefore,  steady  three- 
dimensional  numerical  simulations  of  this  advanced  combustor  concept  are  performed  following 
a  multiphase,  Reynolds-averaged  Navier-Stokes  (RANS),  and  flamelet/progress  variable  (FPV) 
approach  using  FLUENT.  The  specific  objectives  of  this  investigation  are  (1)  to  examine  the 
non-reacting  flow  field  in  terms  of  streamlines  and  g-loads,  (2)  to  inspect  the  reacting  flow  field 
in  terms  of  streamlines,  temperature,  g-loads,  and  equivalence  ratio  (<J>),  (3)  to  determine  the 
effect  of  averaged  centrifugal  force  on  flame  length  and  g-loads,  (4)  to  indirectly  examine  the 
effect  of  centrifugal  force  on  the  turbulent  flame  speed,  and  (5)  to  quantify  the  effect  of 
centrifugal  force  on  exit  temperature  profiles.  Although  there  are  many  more  quantities  of 
interests  not  all  can  be  addressed  in  this  research. 

Physical-Numerical  Procedure 

The  commercial  code  FLUENT  [1]  has  been  utilized  and  has  been  customized  for  the 
simulations  presented  here.  The  fuel  is  injected  as  a  liquid  spray  through  an  injector  and  discrete 
phase  is  used  to  model  the  transition  of  liquid  fuel  spray  into  gas-phase  employing  both  Eulerian 
and  Lagrangian  frames  of  reference.  Details  of  the  non-adiabatic,  non-premixed  FPV  model  are 
presented  here  such  as  gas-phase  formulation,  in-house  effusion  cooling  sources,  geometries  and 
boundary  conditions. 
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2.1.2  Gas-Phase  Numerical  Model 


The  steady  three-dimensional  governing  equations  of  continuity,  momentum,  turbulence,  total 
enthalpy,  mixture  fraction,  mixture  fraction  variance,  and  un-normalized  reaction  progress 
variable  are  solved  using  the  coupled  pressure-based  solver  of  FLUENT  [12].  Turbulence  is 
modeled  using  the  Realizable  k-s  RANS  governing  equations.  Turbulent  kinetic  energy  is  only 
generated  due  to  mean  velocity  gradients.  Scalable-wall  functions  [1,12]  are  used  to  determine 
the  near  wall  flow  velocity.  The  un-normalized  progress  variable  variance  is  obtained  through  an 
algebraic  equation  [12].  The  governing  equations  are  discretized  using  second-order  upwind 
scheme  [1].  The  gradients  and  derivatives  of  the  governing  equations  are  computed  using  the 
Least-Square  Cell-Based  Gradient  [1],  which  is  second-order  spatially  accurate.  The  discrete 
phase  was  set  to  interact  with  the  continuous  phase  every  20  iterations.  The  continuity, 
momentum,  total  enthalpy,  and  the  mixture  fraction  equations  contain  source  terms  that  couple 
the  continuum  phase  with  the  discrete  phase.  The  Eulerian  governing  equations  in  differential 
notation  can  be  represented  in  general  form  as  follows, 


d(£U0 )  d(pv0)  3(pw0)  _  d_  /  0  d0\  d_  /  d<P\  d_  /  0  d<P\  0 

dx  dy  dz  ~  dx  V  dx)  +  dyV  dy)  +  dz  V  dz )  + 


(1) 


Depending  on  the  value  of  <l>,  this  equation  represents  the  continuity,  momentum,  turbulent 
kinetic  energy,  eddy  dissipation  rate,  total  enthalpy,  mixture  fraction,  mixture  fraction  variance, 
and  un-normalized  reaction  progress  variable  equations,  as  shown  in  Table  1.  This  table  also 
describes  the  transport  coefficient  T4*  and  the  source  term  for  each  equation.  Further  details 
regarding  the  model  constants,  and  p-shape  PDF  can  be  found  in  Refs.  [1-3].  Source  and  sink 
terms  in  the  continuity  and  enthalpy  equations  are  used  to  model  effusion  cooling.  These  are 
indicated  by  the  first  term  in  the  So  column  of  Table  1 . 

Discrete  Phase  Model 


Note  that  the  second  terms  in  the  So  column,  the  last  terms  in  the  momentum  equations  and  the 
only  term  in  the  mixture  fraction  equations  are  the  sources  for  the  DPM  interaction  with  the 
Eulerian  phase.  Briefly,  the  fuel  liquid  droplet  trajectory  is  calculated  by  integrating  the  force 
balance  on  the  droplet  by  equating  the  particle  inertia  with  the  drag  force  and  considering 
turbulent  dispersion.  The  droplets  evaporate  when  the  droplet  temperature  (Tp)  is  higher  than  the 
saturation  temperature  (Tiaf)  and  lower  than  the  boiling  point  (Tbp)  of  JP8.  Under  evaporation  the 
JP8  liquid  particle  receives  heat  from  the  continuum  phase  while  transferring  mass  to  the  gaseous 
phase.  For  more  details  refer  to  Refs.  [3,  16,  17]. 


6 

Approved  for  public  release;  distribution  unlimited. 


Table  1.  Variables,  Transport  Coefficients,  and  Source  Terms  Appearing  in  the  Governing  Equation 
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2.1.3  Flamelet  Chemistry 

N-dodecane  is  used  as  a  single  surrogate  component  for  JP8.  The  flamelets  were  calculated  by 
solving  the  laminar  n-C^fWair  counterflow  non-premixed  flame  equations  in  a  mixture  fraction 
space  [1,2]  using  the  JetSurf-1.0  [1]  chemical  reaction  mechanism  containing  194  species  and 
1459  Arrhenius  reactions.  Several  flamelets  for  various  scalar  dissipation,  x,  rates  were 
computed. 
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2.1.4  Un-Normalized  Progress  Variable,  Yc 


The  dimensionless  C-progress  variable  previously  used  [1-3]  is  replaced  by  the  un-normalized 
reaction  progress  variable,  Yc.  The  relationship  between  the  C-progress  variable  and  the  Yc  is  as 
follows. 


C  = 


^eo  +  ^c02 

C(/)+C2(/) 


(2) 


Yc  is  defined  as  the  summation  of  the  CO  and  CO2  mass  fractions.  Note  that  this  definition  is  not 
unique  and  could  have  also  been  defined  in  many  other  ways.  The  C  progress  variable  is  the  ratio 
of  local  Yc  and  the  local  equilibrium  Yc.  It  is  more  advantageous  to  solve  for  Yc  than  for  C.  This 
is  because  the  latter  requires  modeling  of  the  source  terms  based  on  the  laminar  flame  speed  and 
Zimont  model.  On  the  contrary,  the  former  transport  equation  can  achieve  closure  of  the  source 
terms  by  recalling  the  CO  and  CO2  net  reaction  rates  from  the  FGM.  This,  in  turn,  provides 
another  benefit.  That  is,  the  source  term  in  Yc  is  not  just  now  dependent  on  C  but  also  on  f.  This 
implies  that  when  C  =  1,  additional  air  dilution  can  actually  quench  the  flame  (i.e.,  Yc  will  be 
reduced  and  C  will  fall  below  unity). 

2.1.5  Turbulence-Chemistry  Interaction 

The  turbulence-chemistry  interaction  is  modeled  using  probability  density  functions  (PDFs), 
which  can  be  thought  as  the  amount  of  time  a  fluid  particle  spends  in  the  vicinity  of  a  state, 
determined  by  mixture  fraction  (f),  progress  variable  (c),  and  total  enthalpy  (H).  The  FGM’s 
variable  <t»  takes  the  values  of  p,  Yh  T,  and  SFR.  The  triple  integral  of  the  product  of  <t>  and  its 
PDFs  gives  the  averaged  value  of  <t>  in  the  three-dimensional  spatial  domain.  This  general 
equation  is  presented  below. 


f,1  <Kf,  c,  H)PDFa  c,  H)dfdcdH  (3) 

Note  that  f  and  c  are  integrated  from  0  to  1,  whereas  H  is  integrated  from  the  user’s  input  of  Hgain 
to  Hioss-  The  latter  is  needed  for  non-adiabatic  cases  like  the  current  model.  By  decomposing  the 
PDF  into  three  independent  PDFs  the  following  equation  is  obtained. 

$  =  f*  f*  <D(f,  c,  H  )PDF(f)PDF(c)PDF(H)dfdcdH  (4) 

For  the  mixture  fraction  (f)  and  progress  variable  (c),  P  presumed  PDFs  are  used.  These  requires 
second  moments  (i.e.,  variances)  for  f  and  c.  The  mixture  fraction  variance  (f  2)  is  obtained  by  a 
solving  transport  equation  as  shown  in  Table  1.  The  c-progress  variable  variance  is  obtained 
from  the  following  algebraic  consideration. 


v  v-L—Kyic) 

C  =  7^72  =  7 - ^ - rz  (5) 

ft  )  (C(/)+C2W) 

For  the  PDF(H),  a  5  presumed  probability  function  is  used.  This  turns  the  above  equation  into  the 
following. 
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d>  =  £  £  4>(f,  c,  H,)P(f)P(c)S(H  -  H)dfdc 


(6) 


Now  the  thermochemical  table  can  be  constructed  for  the  variables  <D  =  p,  Yj,  T,  and  SFR  as 
function  of  f,  c  and  H.  The  multiple  flamelets  scalar  dissipation  rates,  %,  are  dropped  out  by 
calculating  c  along  a  flame.  Note  below  that  Yj  is  a  function  of  only  f  and  c  within  the  flame 
region  (0  <  c  <  1  )  in  order  to  reduce  the  number  of  flamelets  computed.  However,  Yj  is  also  a 
function  of  H  at  equilibrium  conditions  (c  =  1  ). 


=  /01/01p(f,c,H)P(f)P(c)6(H  -  H)dfdc 

(7) 

=  /0L/0Ln(f,C)P(f)P(C)dfdc,0  <  c  <  1 

(8) 

=  /01/0Vi(f,H)P(f)5(H  -  H)dfdc,c  =  1 

(9) 

=  fl  £  T(f,  c,  H  )P(f)p(c)5(H  -  H)dfdc 

(10) 

$fr  =  Iq  Iq  SFR(f,  c)(3(f)P(c)dfdc 

(11) 

Note  that  SFR  (f,  c)  is  just  the  summation  of  the  net  reaction  rates  of  CO  and  CO2  (i.e.,  SFR  = 

6)co  +  d>co  ).  Moreover,  although  the  transport  equations  solve  for  Yc  this  value  has  to  be 
converted  to  C  through  Eq.  2  in  order  to  interact  with  the  thermochemical  table.  When  local 
strain  rate  weakens  the  turbulent  flame  in  the  flow  field  the  value  of  C  is  less  than  unity  for  a 
given  f.  If  local  strain  rate  is  excessive  the  C  is  zero  and  the  flame  quenches  locally. 

2.1.6  Thermal  and  Transport  Properties 

The  default  values  for  mixture  molecular  viscosity  (p)  and  thermal  conductivity  (X)  in  FLUENT 
are  constant.  This,  however,  does  not  correspond  to  the  reality  of  the  current  conditions. 
Combustion  converts  the  reactant’s  enthalpy  of  formation  into  sensible  enthalpy,  raising  the  flow 
temperature.  This,  in  turn,  decreases  the  density  (p)  up  to  six  times,  while  accelerating  the  flow 
velocity  through  the  flame  and  increasing  p  and  X  The  overall  effect  of  combustion  is  to  locally 
lower  the  Reynolds  number.  Hence,  the  molecular  transport  properties  become  relatively 
important  as  their  turbulent  counterparts  in  the  vicinity  of  the  flame-front  and  high  temperature 
regions.  Therefore,  we  developed,  debugged,  compiled  and  linked  a  C-programming  user- 
defined  function  (UDF)  to  account  for  the  mixture  temperature-  and  species-dependent  p  and  X. 
The  code  reads  the  transport  properties  for  all  species  involved  in  the  flamelet  generation.  Then, 
it  compares  these  species  with  those  available  in  the  PDF  look-up  table  (i.e.,  20  most  abundant 
species  stored  for  current  results).  The  code  automatically  extracts  the  species  molecule  index,  L- 
J  potential,  L-J  collision  diameter,  polarizability  and  the  rotational  relaxation  number  [23].  At  the 
beginning  of  each  iteration,  the  code  computes  the  individual  species  p,  and  X  using  the  standard 
kinetic  theory  [1,  2]  and  the  kinetic  theory  approach  given  by  Wamatz  [22, 26],  respectively.  The 
mixture  molecular  transport  properties  needed  in  the  equations  in  Table  1  are  obtained  from  p, 
and  X  using  Wilke  semi-empirical  equations  [22,  27].  The  in-house  UDF  was  compared  with 
calculations  performed  using  NASA  CEA  [28]  and  they  compared  remarkably  well  when 
comparing  p  and  X  as  a  function  of  §  (not  shown).  This  comparison  also  revealed  that  if 
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temperature-dependent  air  properties  were  to  be  used  alone,  instead  the  error  in  X  would  be  up 
50%  at  fuel-rich  conditions.  Furthermore,  the  transport  kinetic  theory  formulation  is 
computationally  expensive,  slowing  down  iterations  about  5  times.  Thereby,  to  mitigate  this 
drawback  we  blended  the  transport  kinetic  theory  with  air  temperature-dependent  properties.  At 
low  temperatures,  below  850K,  the  latter  transport  properties  would  be  used  whereas  at  higher 
temperature  the  former  transport  properties  would  be  utilized.  This  only  slowed  down  the 
calculations  about  2.5  times.  Finally,  the  mixture  specific  heat  capacity  (cp)  essential  for 
calculating  the  temperature  from  the  enthalpy  equation  is  obtained  from  a  mixing  law  (i.e., 

Cp  =  Yii^i  '  Cp,i)- 

2.1.7  Boundary  Conditions 

Figure  1  presents  the  computational  domain  and  boundary  conditions.  The  UCC-HGC  is 
modeled  as  a  60°  rotational  periodic  sector.  The  geometry  is  represented  by  pressure  inlet,  mass 
flow  inlet,  and  wall  boundary  conditions.  The  UCC-HGC  is  fed  by  a  plenum  associated  with  the 
mainstream  flow  and  the  CC  driver  jets.  Therefore,  this  plenum  is  modeled  as  pressure  inlet 
boundary  conditions.  The  comer  cooling  jets  are  also  modeled  as  pressure  inlet  boundary 
condition.  The  outflow  is  modeled  with  mass  flow  outlet  boundary  condition.  The  reason  for 
using  inlet  pressure  and  mass  outlet  boundary  conditions  is  because  the  various  UCC-HGC 
configurations  studied  here  are  compared  at  constant  mass  flow  rate,  operating  pressure,  and 
inlet  temperature,  while  the  mass  flow  split  between  mainstream  and  CC  driver  jets  is  varied  at 
constant  global  equivalence  ratio  (<(>Giobai)-  The  rest  of  the  walls  are  modeled  as  impermeable 
adiabatic  wall  boundary  conditions.  The  operating  pressure  of  the  combustor  is  1,945  kPa.  The 
fuel  is  injected  as  liquid  from  a  cone  injector  located  at  the  centroid  of  the  outer  surface  of  the 
cavity-in-cavity  (CIC)  feature.  The  injector  is  normal  to  the  CIC  and  CC  and  its  half  angle  of 
injection  is  set  to  0i/2  =  45°  with  a  logarithmic  Rosen-Rammler  droplet  distribution  with  spread 
factor  of  3.5  varying  from  0.1  to  100  pm.  The  average  diameter  is  20  pm.  The  fuel  droplets  can 
escape  the  inlet  and  outlet  boundaries,  whereas  they  bounce  off  the  walls.  The  boundary 
conditions  for  the  continuum  and  discrete  phase  governing  equations  presented  before  are  as 
follows: 


Pressure  Inlet  Boundary  Condition 

P  =  134  kPa  (gauge),  T  =  810.927  K,YC  =  0,f  =  f~2  =  0,u'  =  0.05u,lturb  =  DH 

Mass  Flow  Outlet  Boundary  Condition 
ks 

m  =  0.564  — (per  sector) 

Wall  Boundary  Condition 

u  =  v  =  w  =  0,  VnYc  =  Vnf  =  Vnf"2  =  VnH  =  0,  k  =  <e  =  0,  vp.nw  =  -vp.  nw,  and  vp.  tw  = 
vp-tw 

Discrete  Phase  Injections 

mp  =  0.0189  kg/s  (per  injector),  |v^|  =  93m/s,Tp  =  296.8  K, 01/2  =  45°, dpmin  = 

0.1  |im,  dp  avg  20  |im,  dp  max  100  [im 
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Figure  1.  Computational  Domain  and  Boundary  Conditions 

(The  types  of  boundary  conditions  are  color  coded  in  red,  representing  pressure  inlet  boundary  conditions,  green, 
representing  mass  flow  outlet,  and  gray,  representing  impermeable  adiabatic  wall  boundary  conditions.  The  orange 
arrows  indicate  the  locations  at  which  the  injector  is  mounted  onto  the  circumferential  cavity  (CC).  The  various 
components  of  the  UCC-HGC  are  labeled.  The  labels  are  colored  with  the  applied  boundary  condition.) 

The  mass  and  enthalpy  sink  terms  that  appear  in  Table  1  are  applied  to  the  lower  surface  of  the 
plenum  adjacent  to  the  CC  jets.  Then,  there  are  equivalent  mass  and  enthalpy  source  terms 
adjacent  to  the  outer  CC  surface  as  well  as  the  outer  and  side  surfaces  of  the  CIC.  Here  the  sink 
and  source  terms  are  balanced.  These  sink  and  source  terms  mimic  the  fluid  flow  from  the 
plenum  volume  through  the  effusion  cooling  holes  of  the  outer  CC  surface  of  the  cavity.  In 
addition,  there  are  mass  and  enthalpy  source  terms  applied  to  the  forward  and  aft  surfaces  of  the 
CC  as  well  as  to  the  outer  liner  immediately  downstream  the  CC  but  ends  near  the  trailing  edge 
of  the  IGV.  These  mass  source  terms  model  the  air  effusion  cooling  through  each  of  these 
surfaces. 

Mesh 

We  performed  global  mesh  sensitivity  analysis  because  there  are  nine  geometries  to  consider. 

Our  approach  was  the  following.  We  created  two  mesh  sizes  for  the  nine  configurations  reported 
here,  viz.,  2.5  and  10  million  cells.  Non-reacting  flow  simulations  were  performed  for  these 
mesh  size  and  the  discharge  coefficient  was  computed  for  the  core  and  CC  jets.  Results  indicated 
that  the  discharge  coefficient  was  nearly  independent  of  mesh  size.  The  discharge  coefficient  for 
the  core  flow  and  the  CC  jets  was  in  the  vicinity  of  0.77  and  0.69,  respectively.  We  decided  to 
use  the  10  million  cell  mesh  to  avoid  any  grid  sensible  results.  The  mesh  is  unstructured,  formed 
by  tetrahedral  elements  in  the  core  and  prism  elements  for  the  boundary  layers.  The  minimum 
and  maximum  cell  sizes  for  the  tetrahedral  elements  are  0.50  and  1.50  mm,  respectively.  Ten 
inflation  layers  with  initial  height  of  0.1  mm  are  used  for  the  prismatic  boundary  layers.  We 
decided  to  use  scalable  wall  functions  to  address  the  shear  layer  near  the  walls.  The  mesh 
maximum  stretch  factor  is  1.2. 
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Figure  2.  FLA  Drawing  of  the  Cavity  Volume 

(The  figure  shows  the  driver  jet  tangential  velocity  (—(usin(a))  +  v(cos(a))J  and  its  angle  of  inclination  (a)  as 
well  as  the  general  representation  for  calculating  the  flow  tangential  velocity  (— ( usin(6 ))  +  v(cos(d)))  in  the 
domain.  The  Cartesian  (t,  j)  and  cylindrical  (r,  0)  unit  vector  are  displayed.) 

2.1.8  Configuration  Matrix 

The  steady  volumetric  centrifugal  force  on  a  fluid  particle  in  the  cavity  flow  is  depicted  in  Figure 
2  and  can  be  summarized  as: 

Fc  =  p  ■  (g  -  loads)  -g  and  g  -  loads  =  (12) 

Note  that  the  numerator  is  just  the  tangential  velocity  everywhere  in  the  domain  defined  by  the 
Cartesian  spanwise  (a)  and  transverse  (v)  velocities  and  by  the  circumferential  coordinate,  0.  The 
two  parameters  that  were  chosen  to  vary  centrifugal  force  were  the  jet  flow  inclination  angle 
(a;et)  and  the  cavity  mass  flow  rate  ( mcc ).  It  is  expected  that  the  g-load  should  increase  with 
decreasing  ajet  (cf.  Eq.  12)  and  with  increasing  rhcc.  The  three  chosen  <Xjet  are  represented  in 
the  configuration  identification  names  as  numbers  and  are  as  follows:  1  is  20.0°,  2  is  32.5°,  and  3 
is  45.0°.  The  cavity  mass  flow  rate  ( rhcc )  is  varied  by  changing  the  diameter  of  the  driver  jet 
holes  because  the  pressure  drop  across  the  holes  is  constant  for  all  cases.  In  order  to  maintain 
constant  total  air  flow  and  global  equivalence  ratio  (^Global),  the  mainstream  mass  flow  must 
decrease  with  increasing  mcc.  Thus,  a  blockage  plate  with  vertical  slots  as  shown  in  Figure  1  is 
used  to  meter  this  mainstream  flow.  By  reducing  the  core  flow  with  higher  metering  plate 
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blockage,  we  increase  the  diameter  cavity  driver  jets  so  that  the  combination  of  the  core  flow  and 
the  cavity  driver  jet  flow  is  maintained  nearly  constant.  Therefore,  the  resulting  parameter  is  a 
CC-to-core  area  ratio  where  the  core  area  is  the  open  area  through  the  blockage  plate  and  the  CC 
area  is  the  total  area  of  the  driver  jets.  The  three  chosen  ratios  are  represented  in  the 
configuration  identification  names  as  letters  and  are  as  follows:  “a”  is  0.22,  “b”  is  0.41,  and  “c” 
is  0.74.  Note  that  increasing  the  CC-to-core  area  ratio  directly  increases  mcc/mC0re,  and 
increases  the  cavity  loading  while  leaning  out  the  CC.  There  are  a  total  of  nine  configurations 
investigated  here,  and  they  are  shown  in  Table  2. 


Table  2.  Configuration  matrix 


CC-to-core  Area  Ratio 

0.22 

0.41 

0.74 

20.0° 

la 

lb 

lc 

32.5° 

2a 

2b 

2c 

45.0° 

3a 

3b 

3c 

Table  Note:  There  are  three  CC  jet  angles  (Ujet)  and  for  each  there  are  three  CC-to-core  area  ratios  for  a  total  of 
nine  configurations.  The  alphanumeric  values  in  the  table  are  the  configuration  identification  names. 

2.1.9  Results  and  Discussion 

In  this  section  we  initially  discuss  the  non-reacting  flow  field.  Subsequently,  we  examine  the 
reacting  flow  field  for  various  UCC-HGCs.  Then,  we  determine  the  effect  of  averaged 
centrifugal  force  on  g-loads  and  flame  length.  Next,  we  inspect  the  effect  of  centrifugal  force  on 
the  turbulent  flame  speed.  Finally,  we  quantify  the  effect  of  averaged  centrifugal  force  on  exit 
temperature  profiles. 

Non-reacting  Flow  Field 

The  non-reacting  flow  field  is  examined  in  terms  of  streamlines  and  g-loads.  Figure  3  and  Figure 
4  present  the  streamlines  and  area-averaged  axial  g-load  profiles,  respectively.  Because  the 
streamline  flow  structure  for  the  various  configurations  under  non-reacting  flow  conditions  is 
qualitatively  similar,  Figure  3  only  illustrates  the  flow  field  for  configuration  la.  Note  here  that  a 
plenum  feeds  the  CC  jets  and  some  (blue)  streamlines  end  at  the  CC  outer  surface  because  of  the 
mass  and  enthalpy  sink  terms  imposed  here  to  mimic  transpiration.  Within  the  CC  the  flow  field 
is  very  stratified  with  the  circumferential  flow  marked  by  blue  streamlines  circulating  tangent  to 
the  outer  cavity  wall.  The  mainstream  flow  (red)  does  not  enter  into  the  CC.  The  cooling  jets  are 
washed  out  and  remain  near  the  liner  outer  surface. 

Figure  4  clearly  illustrates  that  there  are  multiple  g-load  peaks  that  are  consistent  and  that  are 
independent  of  the  configuration.  The  two  rows  of  CC  air  driver  jets  produce  two  g-load 
maxima.  The  magnitude  of  the  pair  of  g-load  peaks  in  the  CC  varies  chiefly  as  function  of 
rhcc/ihcore  and  secondly  as  a  function  of  jet  angle  (ajet).  As  a  increases,  the  CC  g-load 
decreases  from  1283  to  916g’s  for  series  a,  from  2661  to  1809g’s  for  series  b,  and  from  4584  to 
3400g’s  for  series  C.  Finally,  the  effect  of  <Xjet  becomes  more  pronounced  at  higher  rhcc/rhC0re. 
The  trends  seen  in  the  non-reacting  flow  results  are  consistent  with  the  trends  inferred  by  Eqs.  12 
and  Figure  2. 
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Figure  3.  FLA  View  of  Non-reacting  Flow  Streamlines 

(The  non-reacting  flow  streamlines  are  emanating  from  the  mainstream  plenum  (red),  the  CC  plenum  (blue),  and  the 
cooling  jets  (green)  for  the  configuration  la  discussed  in  the  context  of  Table  2.) 

Reacting  Flow  Field 

Now  we  discuss  the  reacting  flow  field  for  various  UCC-HGC  configurations.  Figure  5  and 
Figure  6  present,  respectively,  the  forward  looking  aft  (FLA)  and  side  views  of  the  reacting  flow 
streamlines  emanating  from  the  mainstream  plenum  (red),  the  CC  plenum  (blue),  and  the  cooling 
jets  (green),  ft  is  clear  from  Figure  5  that  the  presence  of  the  flame  and  the  resulting  volumetric 
expansion  has  somewhat  distorted  the  (strong)  uniform  CC  bulk  swirl  that  was  present  under 
non-reacting  flow  conditions  (cf.  Figure  3  and  Figure  5).  However,  the  plenum  that  feeds  the  CC 
jets  does  not  seem  to  be  affected  by  combustion  and  streamlines  still  end  at  the  CC  outer  wall 
due  to  the  mass  and  enthalpy  sink  terms  that  mimic  effusion  cooling.  All  configurations  exhibit 
nearly  negligible  mainstream  flow  entraining  into  the  CC  like  their  non-reacting  flow 
counterparts  (cf.  Figure  3).  In  addition,  Figure  5  shows  that  as  mcc /mcore  increases  the  cavity 
streamlines  tend  to  initially  stay  close  to  the  cavity  outer  wall  surface  before  turning  radially 
inward  towards  the  core  flow.  This  would  likely  contribute  to  higher  centrifugal  force,  locally, 
along  the  cavity  outer  wall  surface.  The  large  regions  in  the  cavity  without  streamlines  in 
configurations  2b,  and  2c  show  that  large  recirculation  zones  may  exist  near  the  CC  outer  wall. 

The  side  view  of  the  UCC-HGC'  configurations  in  Figure  6  shows  details  of  the  axial  distribution 
of  the  cavity  driver  air.  For  all  configurations,  the  swirling  cavity  flow  tends  to  migrate  to  the 
downstream  (aft)  wall  of  the  cavity  before  exiting  into  the  mainstream.  This  may  be  due  to  the 
dynamic  pressure  exerted  by  the  mainstream  flow  on  the  shear  region  between  the  two  flows.  As 
a  increases  from  configurations  1  to  3,  the  driver  jet  air  penetrates  deeper  into  the  CC. 
Qualitatively,  it  appears  that  with  increasing  a  the  CC  jet  (blue)  streamlines  mixed  better  with 
the  mainstream  and  cooling  flow  as  illustrated  by  Figure  6.  It  appears  that  there  is  no  significant 
trend  of  riicc/rhC0re  on  the  exit  streamlines. 
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Figure  4.  Axial  Area-averaged  g-load  Profiles  for  Various  Non-reacting  Configurations 

(The  configurations  are  discussed  in  the  context  of  Table  2.  The  region  of  interest  is  from  the  CC  forward  wall  to  the 
CC  aft  wall.  The  magenta  arrows  indicate  the  location  of  the  CC  air  driver  jets.) 


15 

Approved  for  public  release;  distribution  unlimited. 


Figure  5.  FLA  View  of  Reacting  Flow  Streamlines 

(The  reacting  flow  streamlines  emanate  from  the  mainstream  plenum  (red),  the  CC  plenum  (blue),  and  the  cooling 

jets  (green)  for  the  configurations  discussed  in  Table  2.) 
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Figure  6.  Side  View  of  Reacting  Flow  Streamlines 

(The  reacting  flow  streamlines  are  emanating  from  the  mainstream  plenum  (red),  the  CC  plenum  (blue),  and  the 
cooling  jets  (green)  for  the  configurations  discussed  in  the  context  of  Figure  5.) 
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Figure  7.  Isometric  View  of  Temperature  Contours  at  the  Centerplane  and  the  Periodic  Boundary 

Condition 

(The  contour  views  are  for  the  configurations  discussed  in  the  context  of  Figure  5.  The  units  are  in  Kelvin.) 

Temperature  contours  at  two  circumferential  locations,  viz.,  centerplane  and  periodic  boundary, 
are  illustrated  in  Figure  7.  Plenum  air  enters  the  flow  domain  at  the  mainstream  and  cavity 
pressure  inlets.  Temperature  rise  begins  within  the  circumferential  cavity,  and  the  temperature 
continues  to  rise  downstream  of  the  CC.  The  temperature  in  the  cavity  is  generally  highest 
toward  the  centerplane.  The  air-cooled  walls  tend  to  keep  hot  products  away  from  the  surfaces. 

It  is  evident  from  this  figure  that  configurations  with  ctjet= 20.0°  (la,  lb,  and  lc)  have  relatively 
higher  temperatures  in  the  CC  than  the  other  configurations.  Interestingly,  configurations  2b  also 
attained  relatively  higher  temperature  region  within  the  CC.  Lastly,  there  is  no  discemable  effect 
of  mcc/mC0re  on  the  temperature  contours. 

The  temperature  contours  presented  in  Figure  7  can  be  better  analyzed  by  conducting  data 
reduction.  This  is  done  by  plotting  the  area- averaged  axial  temperature  profiles  in  Figure  8. 
Because  the  local  equivalence  ratio,  <|),  is  one  of  the  variables  that  determines  the  flame 
temperature,  the  area-averaged  4>  as  a  function  of  dimensionless  axial  distance  (L*)  is  shown  for 
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all  configurations.  These  results  indicate  that  for  all  the  configurations  the  temperature  in  the  CC 
is  lower  (1300-1800K)  than  the  end  of  the  rig  (ca.  2000K).  Importantly,  the  temperature  does  not 
reach  its  maximum  temperature  by  the  trailing  edge  of  the  vane;  indeed,  it  continues  increasing 
slightly  through  the  expansion  region  over  the  centerbody.  This  indicates  further  burning  is 
occurring  in  this  expansion  region.  The  implication  of  this  result  is  that  a  UCC-HGC  would 
actually  end  at  the  IGV  trailing  edge  and  further  burning  would  continue  towards  the  turbine. 
This  result  would  be  detrimental  for  proper  engine  operation.  The  relatively  lower  area-averaged 
temperature  region  inside  the  CC  is  partly  due  to  the  higher  <j>,  well  above  stoichiometry.  For  all 
configurations  the  temperature  at  the  exit  of  the  rig  is  ca.  2000  K  because  the  (|)Giobai=0.53  for  all 
configurations,  as  indicated  at  L*=l.  This  is  consistent  with  equilibrium  gas  temperature 
calculations  for  Jet-A  type  fuel  at  the  rig  operating  pressure,  inlet  air  temperature,  and  (|)Giobai- 


Figure  8.  Axial  Area-averaged  Temperature  and  <|>  Profiles  as  a  Function  of  Dimensionless  Axial 

Position,  L*. 


Effect  of  Centrifugal  Force  on  Flame  Length 

Of  great  importance  for  the  UCC-HGC  design  is  to  address  the  effect  of  centrifugal  force  on 
flame  length.  Here  the  flame  surface  is  defined  as  the  stoichiometric  iso-surface  holding 
combustion  as  shown  in  Figure  9.  The  iso-surfaces  are  colored  by  temperature.  The  results 
indicate  that  the  flame  extends  from  the  CC  driver  jets  towards  the  trailing  edge  of  the  IGV.  It  is 
important  to  note  that  even  though  the  iso-surface  is  at  constant  equivalence  ratio  the  temperature 
varies  along  the  flame  surface.  Most  flames  are  typically  cooler  at  the  CC,  which  is  consistent 
with  Figure  7  and  Figure  8.  Further  downstream  the  temperature  increases  and  it  is  maximum  at 
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the  tail  of  the  flame  front.  This  occurs  because  by  varying  the  CC  jet  angle  (ajet)  and/or  mcc/ 
iV-core  the  stretch  rate  or  scalar  dissipation  rate  changes.  These  flame  quenching  effects  are  more 
pronounced  in  the  CC  than  further  downstream  the  CC.  Hence,  the  CC  flame  temperature  is 
lower  than  its  corresponding  equilibrium  flame  temperature. 

Qualitatively,  increasing  the  surface  area  and  temperature  of  the  flame  that  can  be  sustained 
within  the  CC  leads  to  lower  flame  length.  The  converse  is  also  true.  Configuration  lc  and  3b 
exhibit  the  shortest  flame  lengths.  They  achieve  this  by  increasing  their  volume  within  the  CC. 
Therefore,  longer  flames  typically  contained  small  flame  volume  within  the  CC.  In  terms  of 
flame  length  as  a  function  of  a  and  mcc /rhcore  no  definite  pattern  is  observed. 


Temperature 


Figure  9.  Side  View  of  the  Stoichiometric  Iso-surface 

(The  iso-surface  is  clipped  by  C  >  0. 5  and  colored  by  temperature  (in  Kelvin).  The  upstream  portion  of  the  UCC- 

HGC  was  truncated  for  convenience  in  this  image.) 

Effect  of  Centrifugal  Force  on  g-loads  and  Turbulent  Flame  Speed 

Since  the  design  of  this  new  rig  that  can  operate  at  higher  Mach  number,  higher  operating 
pressures  and  higher  inlet  temperatures,  while  sacrificing  centrifugal  force,  is  based  on  non¬ 
reacting  flow  consideration  it  is  vital  to  understand  what  the  effect  of  combustion  on  the  g-loads 
is.  Figure  10  displays  the  area-averaged  axial  g-loads  profiles  for  all  reacting  flow  configurations 
under  investigation.  Important  observations  are  as  follows.  First,  the  large  g-load  peaks  at  the  CC 
are  diminished  at  low  Thcc/mC0re  while  enhanced  at  high  mcc/mC0re  (cf.  Figure  4  and  Figure 
10).  Second,  at  intermediate  CC-to-core  air  mass  flow  rates  the  g-loads  are  boosted  at  low  a, 
while  g-loads  are  weakened  at  high  ajet.  Therefore,  there  is  threshold  dependent  on  rhcc/mcore 
and  a.jet  at  which  CC  g-loads  are  boosted  or  weakened.  The  g-loads  are  diminished  because  the 
velocity  vector  deflection  occurs  through  the  flame  front.  This,  in  turn,  enhances  the  radial  and 
axial  velocity  through  thermal  expansion,  while  reducing  the  magnitude  of  the  tangential 
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-loads 


velocity.  On  the  contrary,  when  g-loads  are  enhanced  the  flame  front  is  nearly  perpendicular  to 
the  bulk  swirl  flow  and  the  flame  accelerates  the  velocity  vectors  primarily  in  the  tangential 
direction.  Third,  configuration  lc  reaches  the  maximum  g-load,  ca.  6000g,  which  is  above  its 
non-reacting  value  of  4500g.  Similarly  to  non-reacting  flow  results  which  are  very  consistent 
with  variations  in  CC  jet  angle  (ajet)  and  mcc/mC0re,  the  reacting  flow  results  also  very 
consistent.  However,  reacting  flow  results  suggests  that  ajet=20.0o  and  a;et=45.0°  are  more 
sensitive  to  changes  in  mcc/mC0re  than  <Xjet= 32.5°  is.  In  addition,  high  mcc/mC0re  is  more 
sensitive  to  changes  in  a  than  lower  mcc /rhcore  are. 


Figure  10.  Area-averaged  Axial  g-loads  Profiles  as  a  Function  of  Dimensionless  Axial  Position,  L* 

(The  region  of  interest  is  from  the  CC  forward  wall  to  the  CC  aft  wall.  The  magenta  arrows  indicate  the  location  of 

the  CC  air  driver  jets.) 

Turbulent  flame  speed  has  always  been  a  major  interest  as  discussed  in  the  Introduction  section. 
However,  it  is  not  easy  to  calculate  it  from  very  complicated  geometry  and  turbulent  partially 
premixed  flames.  Therefore,  we  indirectly  examine  the  turbulent  flame  speed  by  inspecting  the 
un-normalized  reaction  progress  variable  source  term.  If  turbulent  flame  speed  increases,  then, 
the  reaction  progress  variable,  Sc,  should  also  increase.  Figure  1 1  illustrates  the  area-averaged 
axial  Sc  for  all  the  configurations  under  investigation.  Evidently,  there  are  two  regions  of  Sc- 
Low  Sc  occurs  within  the  CC  while  very  large  Sc  occurs  immediately  downstream  the  CC.  This 
indicates  that  most  of  the  reaction  occurs  outside  the  CC,  which  is  consistent  with  Figure  7  and 
Figure  9.  By  comparing  Figure  10  and  Figure  1 1,  it  is  suggested  that  the  g-loads  are  enhancing 
Sc  and  potentially,  turbulent  flame  speed,  especially  at  high  rhcc/rhC0re.  This  enhancement 
occurs  through  mean  velocity  gradients,  which  generate  turbulent  kinetic  energy,  that  wrinkle  or 
corrugate  the  flame  surface  area.  This  does  not,  however,  necessarily  mean  that  the  turbulent 
flame  speed  is  various  factors  greater  than  its  corresponding  turbulent  flame  speed  given  the 
same  flow  conditions.  Nevertheless,  Sc  is  lower  in  the  CC  than  in  the  IGV  passage  plausibly 
because  the  stretch  rates  are  high  in  the  CC  as  indicated  by  high  g-loads.  Downstream  of  the  CC 
the  stretch  rates  are  relatively  low  and  Sc  increases. 
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Figure  11.  Area-averaged  Axial  Sc  Profiles  as  a  Function  of  Dimensionless  Axial  Position,  L*. 


Effect  of  Centrifugal  Force  on  Exit  Temperature  Profile 


Finally,  attainment  of  the  desired  temperature  profile  is  of  paramount  importance  because  an 
adverse  temperature  profile  can  damage  the  turbine  blades  and  deteriorate  turbine  performance. 
Since  stresses  are  highest  at  the  turbine  hub  and  seal  materials  need  to  be  protected  at  the  turbine 
blade  tip,  the  temperature  profile  needs  to  be  center-peaked  for  modem  high-performance 
engines  [29],  Figure  12  shows  that  area-averaged  radial  temperature  (Tave)  profiles  just 
downstream  of  the  IGV  trailing  edge.  For  series  a,  the  best  configuration  is  la,  which  has  a  CC 
jet  angle  (ajet)  of  20.0°.  The  exit  temperature  profile  is  nearly  center-peaked.  As  the  mcc/mC0re 
or  ajet  increases  the  center-peaked  profile  tends  to  deteriorate.  The  exit  temperature  profile  is 
less  sensitive  to  a.jet  at  high  rhcc/rhcore  than  it  is  under  low  conditions.  Similarly,  exit 
temperature  profiles  are  more  sensitive  to  rhcc/rhC0re at  small  ajet  than  at  large  ctjet. 


Figure  12.  Area-averaged  Radial  Temperature  (TaVe)  Profiles  just  Downstream  of  the  IGV  Trailing 

Edge 

(The  profiles  are  for  the  configurations  discussed  in  the  context  of  Figure  5.  The  radial  coordinate  is  dimensionless, 

Re.it*-) 
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2.1.10  Conclusions 


A  numerical  investigation  of  reacting  flows  in  an  advanced  high-g  cavity,  Ultra-Compact 
Combustor  (UCC)  concept  is  conducted.  This  UCC  concept  combines  the  primary  and  the 
dilution  zones  of  a  conventional  combustor  with  high  pressure  turbine  inlet  guide  vanes  (IGVs) 
to  shorten  the  length  of  the  combustor  compared  to  a  conventional  design  with  separate 
combustor  and  IGV  components.  The  high-g  cavity  UCC  (UCC-HGC)  design  uses  high  swirl  in 
a  circumferential  cavity  (CC)  wrapped  around  a  main  stream  annular  flow.  The  high  swirl  is 
generated  through  angled  CC  driver  jets,  which,  in  turn,  is  based  on  the  driver  jet  injection  angle 
(a jet)  and  the  pressure  drop  across  the  driver  jets.  This  centrifugal  force  is  varied  by  changing 
the  mcc/mC0re  and  inclination  angle  (ajet)  relative  to  the  cavity  ring  surface.  Steady,  rotational 
periodic,  three-dimensional  numerical  simulations  of  this  advanced  combustor  concept  are 
performed  following  a  multiphase,  Reynolds-averaged  Navier-Stokes  (RANS),  and  nonpremixed 
flamelet  (PVF)  approach  using  FLUENT.  Important  conclusions  are  as  follows. 

1.  Under  non-reacting  flow  conditions  the  results  for  all  configurations  analyzed  here  are  similar.  The  plenum 
feeds  both  the  CC  and  the  mainstream  (core)  flow.  The  driver  jets  imposed  a  very  clear  bulk  swirl  flow  within 
the  CC  and  the  mainstream  flow  does  not  entrain  into  the  CC.  However,  under  reacting  flow  conditions  the 
flame  interacts  with  the  flow  and  the  bulk  swirl  flow  is  no  longer  very  strong,  especially  at  low  CC-to-core  air 
mass  flow  rates  (jhcc/mC0re)  and  larger  ajet.  This  is  plausibly  because  the  flame  deflects  the  incoming  driver 
jet  flow,  enhancing  radial  and  axial  velocity  components  through  thermal  expansion  while  diminishing  the 
tangential  component. 

2.  Under  non-reacting  flow  conditions,  the  maximum  g-load  clearly  decreases  with  increasing  CC  driver  jet  angle 
(i tZjet ).  This  effect  becomes  even  more  pronounced  with  increasing  mcc /mcore .  Nonetheless,  the  g-loads  are 
more  sensitive  to  the  increase  in  rhcc/mcore  than  to  the  variation  in  a.  On  the  other  hand,  under  reacting  flow 
conditions,  the  high  CC  g-loads  are  diminished  at  low  rhcc/mC0re  and  large  ajet.  Again,  this  is  plausibly 
because  the  flame  deflects  the  incoming  driver  jet  flow.  At  high  rhcc/mC0re  and  small  a  the  g-loads  are  instead 
further  enhanced  beyond  those  exhibited  under  non-reacting  flow  conditions.  This  is  possibly  due  that  the  flame 
is  no  longer  deflecting  the  flow  but  instead  it  is  accelerating  in  the  tangential  velocity  because  the  flame  is  now 
more  perpendicular  to  the  bulk  swirl  flow.  This  is  exciting  because  the  g-loads  do  not  necessarily  decrease  with 
reducing  combustor  size  if  large  CC  flow  rates  are  applied. 

3.  Most  flames  are  typically  cooler  inside  the  CC  but  further  downstream  the  temperature  increases  and  peaks  at 
the  flame  tail.  This  suggests  that  flame  quenching  effects  are  more  pronounced  in  the  CC  than  in  the 
mainstream.  Qualitatively,  increasing  the  surface  area  and  temperature  of  the  flame  that  can  be  sustained  within 
the  CC  leads  to  lower  flame  length.  The  converse  is  also  true.  There  appears  to  be  negligible  effect  ofmcc/ 
rhcore  on  flame  length. 

4.  Even  though  g-loads  appear  to  enhance  reaction  progress  variable  source  (Sc)  and,  consequently,  turbulent 
flame  speed,  this  does  not  necessarily  mean  that  the  turbulent  flame  speed  under  g-loads  exceeds  its 
corresponding  turbulent  flame  speed  under  0  g’s. 

5.  In  terms  of  exit  temperature  profiles,  the  configuration  that  exhibits  the  best  center-peaked  profile  is  that  at 
smallest  ajet  and  low  mcc  /mcore .  As  the  CC-to-core  mass  flow  ratio  or  ajet  increases  the  center-peaked 
profile  tends  to  deteriorate.  The  exit  temperature  profile  is  less  sensitive  to  ajet  at  high  mcc /mcore  than  it  is 
under  low  conditions.  Similarly,  exit  temperature  profiles  are  more  sensitive  to  m cc/rhC0re  at  small  ajet  than  at 
large  ajet . 
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2.2  Influence  of  Fuel  Characteristics  on  the  Prediction  of  Lean  Blowout  in  Bluff- 
body  Stabilized  Flames 

2.2.1  Introduction 

Gas  turbine  engines,  both  in  ground-based  industrial  settings  and  in  aviation  settings,  utilize 
bluff-body  stabilized  flames,  often  running  at  lean  fuel-air  ratios  to  extend  the  lifetime  of  the 
equipment  or  to  comply  with  emissions  regulations.  However,  running  the  equipment  at  a  lean 
condition  also  puts  the  system  at  risk  for  lean  blowout,  which  can  result  in  facility  inefficiencies, 
hardware  damage,  and  a  catastrophic  reduction  in  aircraft  performance.  For  this  reason 
understanding  lean  blowout  and  predicting  its  occurrence  is  critical  for  the  operability  of  these 
systems. 

Some  of  the  first  experimental  studies  of  lean  blowout  behavior  in  bluff-body  stabilized  flames 
were  conducted  by  Williams  [1]  in  the  late  1940s  and  DeZubay  [2]  in  the  1950s.  DeZubay’s 
work,  stemming  from  experiments  conducted  on  disk-shaped  flameholders,  focused  on 
correlating  the  fuel-to-air  ratio  at  lean  blowout  with  a  parameter  based  on  the  inlet  pressure  and 
velocity  and  the  characteristic  diameter  of  the  flameholder.  In  his  original  work  in  1950, 
DeZubay  did  not  vary  the  inlet  temperature  of  the  rig;  therefore,  a  modified  DeZubay  correlation 
that  includes  temperature  dependence  has  become  more  widely  used,  with  T  representing  the 
temperature  parameter: 


U*  104 

p  0.95  ^-^0.85  rp\2 


(13) 


In  this  correlation,  velocity  (U)  is  in  ft/s,  pressure  (P)  in  psi,  diameter  (D)  in  inches,  and 
temperature  (T)  in  Rankine.  Although  DeZubay’s  experiment  included  both  rich  and  lean 
blowout,  the  correlation  is  intended  only  for  the  lean  region  of  the  stability  curve.  Figure  1  shows 
the  DeZubay  curve  as  constructed  from  the  data  used  for  the  present  work. 

While  the  DeZubay  curve  is  a  well-known  and  popular  correlation,  Figure  1  demonstrates  some 
of  its  drawbacks.  At  a  given  DeZubay  number,  the  equivalence  ratio  at  lean  blowout  can  vary 
greatly.  For  example,  at  a  DeZubay  number  near  6000,  the  equivalence  ratio  at  lean  blowout 
ranges  from  0.28  up  to  0.75.  In  addition,  there  appears  to  be  a  fuel  effect  that  the  correlation 
does  not  capture.  The  propane  data  follow  a  curve  reasonably  well,  as  indicated  by  the  black  line 
in  Figure  1.  However,  data  from  experiments  run  on  JP-4,  methane,  gasoline,  and  jet  fuels  7721, 
6308,  and  6169  follow  different  patterns.  The  AFRL  jet  fuel  data  is  new  data  from  the  current 
study. 
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Figure  13.  Dezubay  Correlation  with  Data  from  Present  Work 

Following  DeZubay’ s  work  the  lean-blowout  behavior  of  bluff-body  stabilized  flames  continued 
to  be  studied  by  numerous  researchers,  including  King  [3],  Zukoski  and  Marble  [4,  5],  Ozawa 
[6],  Ballal  and  Lefebvre  [7],  and  Plee  and  Mellor  [8].  Most  researchers  agreed  that  lean  blowout 
was  governed  by  a  competition  between  the  chemical  and  fluid-mechanic  timescales,  which  can 
be  ratioed  to  yield  what  is  known  as  a  Damkdhler  number.  However,  researchers  have  debated 
which  timescales  are  critical  to  the  lean  blowout  process. 

According  to  Zukoski  and  Marble  [5],  active  combustion  in  the  region  of  the  recirculation  zone 
is  limited  to  its  shear-layer  boundaries,  with  a  portion  of  the  hot  products  being  re-entrained  into 
the  recirculation  zone.  In  a  stable  flame,  the  fresh  gas  flowing  past  the  bluff  body  is  ignited 
before  it  reaches  the  end  of  the  recirculation  zone,  propagating  the  flame.  Zukoski  and  Marble 
theorize  that  when  the  fresh  gas  fails  to  mix  with  the  hot  products  in  the  recirculation  zone  for  an 
adequate  period  of  time,  the  flame  will  blow  off.  They  identified  the  time  that  the  fresh  gas 
spends  in  the  shear  layer  and  the  ignition  delay  time  as  the  critical  timescales  in  the  blowout 
process.  They  also  identified  the  ratio  of  the  free-stream  velocity  and  the  length  of  the 
recirculation  zone  as  being  representative  of  the  time  that  the  fresh  gas  spends  in  contact  with  the 
hot  combustion  products. 

Plee  and  Mellor  [8]  also  developed  a  model  to  explain  lean  blowout  that  was  based  on  the 
various  timescales  involved,  including  fluid-mechanic  and  chemical  times  as  well  as  droplet- 
evaporation  time  and  fuel-injection  mixing  time  in  the  case  of  heterogeneous  mixing.  The  fluid- 
mechanic  timescale  was  taken  to  be  the  turbulent  mixing  time  in  the  shear  layer  between  the 
relatively  cold  incoming  reactants  and  the  hot,  burned  gases  in  the  recirculation  zone.  This 
timescale  was  approximated  as  the  ratio  of  the  flameholder  width,  which  was  also  considered  to 
be  a  measure  of  the  size  of  the  recirculation  zone,  and  the  velocity  of  the  incoming  gas  at  the  lip 
of  the  flameholder.  The  ignition  delay  time  was  taken  to  be  the  critical  chemical  timescale.  The 
model  developed  by  Plee  and  Mellor  successfully  collapsed  the  lean-blowout  data  from  Ballal 
and  Lefebvre  [7]  to  a  straight  line. 

Radhakrishnan  et  al.  [9]  developed  a  correlation  for  the  blowout  velocity  of  premixed,  bluff- 
body  stabilized  turbulent  flames  by  examining  the  relevant  turbulent  and  chemical  timescales, 
using  a  model  of  turbulence  where  small-scale  structures  are  actually  similar  to  tubes  of 
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vorticity.  Radhakrishnan  et  al.  argued  that  the  Taylor  microscale  is  the  scale  across  which  the 
flame  must  propagate  at  the  laminar  flame  speed,  SL,  and  that  for  the  flame  to  remain  stable,  the 
laminar  flame  propagation  must  occur  in  a  time  less  than  the  characteristic  fluid-mechanic  time. 
Based  on  previous  studies,  including  those  of  Plee  and  Mellor  [8],  Radhakrishnan  et  al.  [9]  took 
the  shear-layer  mixing  time  to  be  the  critical  characteristic  turbulent  timescale  rather  than  the 
residence  time  in  the  recirculation  zone. 

Assuming  isotropic  turbulence  and  that  the  turbulence  intensity  is  proportional  to  the  velocity  at 
the  flameholder  lip  (Uhp),  and  the  integral  length  scale  of  a  single  eddy  is  proportional  to  the 
length  of  the  recirculation  zone  (L),  Radhakrishnan  reduced  the  Damkohler  number  to  Eq  14, 
where  v  is  the  kinematic  viscosity,  xc  the  chemical  timescale,  xe  the  eddy  (fluid  mechanic) 
timescale,  and  SL  the  laminar  flame  speed. 


y 

v  l  j 


(sLy 


(14) 


Using  Eq  14  as  well  as  empirical  correlations  to  estimate  the  laminar  flame  speed, 

Radhakrishnan  et  al.  [9]  developed  a  model  for  the  blowout  velocity  of  bluff-body  stabilized 
flames  that  correlated  well  with  experimental  data  obtained  on  lean  propane  flames.  The  lip 
velocity  was  found  to  be  a  better  parameter  for  determining  stability  than  the  gas  inlet  velocity. 

In  addition,  the  recirculation  zone  length  was  found  to  be  a  more  appropriate  length  scale  than 
the  flameholder  size,  although  the  authors  acknowledged  the  difficulties  in  measuring  the 
recirculation  zone  length. 

Kariuki  et  al.  [10]  developed  a  parameter  based  on  the  work  of  Radhakrishnan  et  al.  [9]  and  then 
plotted  this  value  against  the  flame  power  for  swirl-stabilized  flames;  they  found  that  the  scatter 
was  “considerably  smaller”  than  that  usually  observed  in  the  literature. 

Research  conducted  by  Lieuwen  and  colleagues  [11-13]  has  also  been  focused  on  the  mechanism 
of  lean  blowout  and  the  timescales  that  may  be  critical  to  the  process.  Nair  et  al.  [11]  conducted 
laser-sheet  imaging  studies  and  velocity-field  measurements  using  particle  image  velocimetry 
(PIV)  to  determine  the  transient  dynamics  of  bluff-body  stabilized  flames  near  blowout.  They 
found  that  the  first  stage  of  blowout  was  characterized  by  localized  “holes”  in  the  flame  sheet 
where  the  instantaneous  stretch  rate  exceeds  the  extinction  stretch  rate.  As  the  equivalence  ratio 
decreased  from  this  condition,  the  duration  and  scale  of  these  “holes”  increased.  The  authors 
argued  that  as  blowout  was  approached,  the  flame  entered  a  second  stage  of  blowout  where 
large-scale,  violent  flapping  of  the  flame  front  was  noted,  similar  to  von  Karman-type  shedding. 
They  suggested  that  the  possible  mechanism  for  blowout  is  the  introduction  of  a  “slug”  of 
relatively  cool,  unreacted  mixture  into  the  recirculation  zone  by  the  large-scale  dynamics 
observed  in  the  second  stage.  This  “slug”  would  then  fail  to  ignite  the  incoming  mixture,  and  the 
flame  would  blow  out.  Chauduri  et  al.  [14]  noted  a  similar  progression  of  flame  behavior  near 
blowoff,  where  shear-layer  extinction  and  recirculation  zone  burning  are  precursors  to  blowoff 
for  premixed  flames. 

Husain  [12]  attempted  to  develop  a  Damkohler  number  that  could  adequately  model  lean 
blowout  in  a  large  set  of  data  obtained  from  past  and  recent  literature,  noting  that  almost  every 
past  researcher  has  concluded  that  a  competition  between  some  fluid-mechanic  timescale  and  a 
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chemical  timescale  is  critical  in  explaining  the  blowout  process.  Using  CHEMKIN  software, 
Husain  calculated  the  minimum  residence  time  of  a  well-stirred  reactor  for  which  a  majority  of 
the  reactants  are  consumed  (tpsr)  and  used  this  as  the  chemical  timescale  in  the  Damkohler 
number.  For  the  fluid-mechanic  timescale,  Husain  assumed  D/U,  where  U  is  the  lip  velocity, 
and  D  is  the  characteristic  diameter  of  the  bluff  body.  Husain  claims  a  “reasonable  grouping”  of 
the  data  based  on  a  plot  of  his  Damkohler  number  as  a  function  of  Reynolds  number,  and  that  the 
Damkohler  number  shows  variation  with  Reynolds  number.  However,  the  dependence  on 
Reynolds  number  is  clearer  within  a  few  individual  datasets  than  the  data  set  as  a  whole,  and  the 
plot  itself  is  on  a  log  scale. 

While  bluff-body  flames  have  been  studied  extensively  for  decades,  a  consensus  on  the  physical 
and  chemical  processes  that  lead  to  blowout  has  not  been  reached.  The  conclusions  of 
Shanbhogue  et  al.  [13]  imply  that  the  mechanism  relies  on  ignition  processes.  Both 
Radhakrishnan  et  al.  [9]  and  Kariuki  et  al.  [10]  had  some  success  correlating  the  equivalence 
ratio  at  blowout  to  flame  speed,  which  would  suggest  that  blowout  is  not  governed  by  the 
ignition  timescales.  Shanbhogue  et  al.  [13]  also  suggest  that  as  the  flame  nears  blowout,  it 
develops  a  sinuous  instability,  which  entrains  cold  reactants  into  the  flame,  disrupting  the 
ignition  process.  However,  it  is  unclear  whether  all  bluff-body  flames  exhibit  this  sinuous 
structure  before  blowout  or  whether  it  has  a  primary  influence  on  the  blowout.  For  example, 
Khosla  et  al.  [15]  performed  a  large  eddy  simulation  (LES)  analysis  on  bluff  bodies  with  and 
without  small  tabs.  The  tabs  effectively  broke  up  the  von  Karman  vortex  shedding,  and  yet  both 
flameholders  blew  out  at  the  same  premixed  equivalence  ratio.  In  addition,  a  study  by 
Huelskamp  et  al.  [16]  evaluated  the  effect  of  shape  on  blowout  by  correlating  the  cold-flow 
Strouhal  number  of  various  bluff-body  shapes  to  the  equivalence  ratio  at  lean  blowout.  The 
Strouhal  number  did  not  appear  to  have  a  first-order  effect  on  lean  blowout,  indicating  that  the 
shape  and  the  vortex  shedding  behavior  did  not  govern  the  mechanism  of  blowout. 

The  present  work  employed  experimental  data  taken  at  the  Air  Force  Research  Laboratory  and 
that  collected  from  many  of  the  papers  reviewed  above  to  develop  a  correlation  for  predicting 
lean  blowout  using  a  least-squares  curve-fit  method.  The  laminar  flame  speed  and  ignition  delay 
time  were  calculated  for  current  experimental  conditions  as  well  as  for  selected  data  sets  from 
the  literature  using  the  chemical  kinetics  software  Cantera,  and  the  results  were  incorporated  into 
the  correlations. 

This  study  will  serve  two  purposes:  1)  to  provide  an  accurate,  practical  method  of  predicting 
lean  blowout  for  designers  and  modelers  and  2)  to  provide  insight  into  the  critical  parameters  and 
timescales  that  govern  the  blowout  process  by  examining  the  significance  of  each  parameter 
included  and  the  physical  and  chemical  processes  that  it  may  affect. 

2.2.2  Experimental  Setup 

Experiments  were  performed  in  two  combustion  laboratories  at  the  Air  Force  Research 
Laboratory  (AFRL)  at  Wright-Patterson  Air  Force  Base.  Bluff-body  flameholding  tests  were 
initially  performed  in  an  atmospheric-pressure  combustion  facility  using  propane  as  a  fuel.  The 
facility  can  deliver  air  inlet  temperatures  of  up  to  500  K  with  four  electric  heaters  with  a  total 
power  capacity  of  358  kW.  Inlet  temperatures  of  up  to  1000  K  can  be  achieved  with  vitiation. 

Air  is  metered  through  orifice  plates  with  an  uncertainty  within  2%,  and  coriolis  meters  measure 
the  fuel  flow  within  0.35%  of  total  flow. 
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The  available  airflow  to  this  facility  can  provide  lip  Mach  numbers  of  up  to  0.3  and  lip  Reynolds 
numbers  up  to  100,000  in  a  test  section  that  is  15.2  cm  high  by  12.7  cm  wide.  Lip  Reynolds 
number  is  calculated  using  the  flameholder  diameter  and  the  velocity  of  the  incoming  gas  as  it 
passes  over  the  flameholder.  Figure  14  is  a  schematic  of  the  test  rig  with  a  v-gutter  bluff  body 
installed.  The  airflow  is  left  to  right  and  passes  through  a  perforated  plate,  which  provides  more 
uniform  velocity  and  temperature  profiles  at  the  inlet  of  the  test  section. 

The  flameholder  trailing  edge  was  58  cm  from  the  inlet  of  the  test  section.  The  flameholder 
spans  the  full  height  of  the  rig,  and  optical  access  through  quartz  windows  is  available  for  high¬ 
speed  imaging.  Gaseous  propane  is  premixed  with  the  air  ~  150  cm  upstream  of  the  test  section. 
The  fuel  distribution  was  validated  using  acetone-seeded  planar  laser-induced  fluorescence 
(PLIF).  The  fuel-air  mixture  is  ignited  with  an  ethylene  torch  behind  the  bluff  body.  Figure  15 
shows  the  rig  as  it  is  installed  in  the  facility. 


Figure  14.  Atmospheric  Pressure  Combustion  Rig  Configuration,  (A)  Top  View;  (B)  Side  View 


Figure  15.  Atmospheric  Pressure  Combustion  Rig 

Experiments  for  this  study  were  also  performed  in  the  High  Pressure  Combustion  Research 
Facility.  Various  jet  fuels,  including  petroleum-derived  JP-8,  tallow-derived  hydro-treated 
renewal  jet  (HRJ)  fuel,  and  camelina-derived  HRJ  fuel,  were  combusted  in  a  vitiated  stream. 

The  molecular  weight  and  hydrogen-to-carbon  ratios  of  these  fuels  were  measured  by  Princeton 
University  [17]. 
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Table  3.  Correlations  Considered  in  This  Study 


Reference 

Lip 

Velocity 

(m/s) 

Flame- 

holder 

Diameter 

(cm) 

Inlet 

Temperature 

(K) 

Inlet 

Pressure 

(atm) 

Reynolds 

Number 

Fuel  Type 

Huelskamp  et 
al. 

6.40-52.1 

0.95-3.8 

290-574 

0.98 

4800-86,000 

Propane 

Ballal  and 
Lefebvre 

10.3-151 

2.0-12.5 

300-575 

0.20- 

0.99 

12,000-506,000 

Propane 

DeZubay 

26.2-215 

0.64-2.5 

306 

0.20-1.0 

4200-345,000 

Hopane 

King 

122-226 

3.8 

700-1033 

0.35- 

0.85 

14,000-87,000 

JP-4 

Potter  and 

Wong 

24.1-233 

0.95-2.5 

300-306 

0.25- 

0.94 

7,300-260,000 

Propane 

Williams 

8.84-108 

0.25-1.3 

300 

1.0 

1,400-89,000 

Propane 

Zukoski 

24.7-368 

0.025-0.64 

339-478 

1.0 

200-57,000 

Methane 

Gasoline 

Chaudhuri  et  al. 

20.1-39 

0.95 

285-295 

1.0 

13,000-25,000 

Propane 

Yamaguchi 

9.14-38.1 

2.4 

293 

1.0 

14,000-60,000 

Propane 

Barrere 

11.9-47.9 

0.50 

290 

1.0 

4,000-16,000 

Propane 

AFRL  (current) 

72.8-149 

3.8 

732-920 

0.68 

20,000-53,000 

Jet  fuels 

The  test  section  is  15.2  cm  wide  and  22.9  cm  tall  with  a  3.8-cm-wide  v-gutter,  and  Mach 
numbers  up  to  0.3  can  be  achieved.  In  this  facility,  the  pressure  in  the  test  section  can  be  varied 
from  0.34  to  1.4  atm.  Heaters  provide  inlet  air  temperatures  of  up  to  670  K,  and  vitiation  can 
raise  the  test  section  inlet  temperature  to  1000  K.  Oxygen  levels  of  the  gas  entering  the  test 
section  can  be  varied  from  17  to  21%,  depending  on  the  chosen  level  of  vitiation;  for  this  study, 
the  points  chosen  were  taken  at  an  oxygen  level  of  ~18%.  The  test  section  fuel  is  injected  ~51 
cm  upstream  of  the  flameholder  through  atomizing  nozzles,  providing  a  premixed,  pre-vaporized 
fuel  condition.  The  fuel  distribution,  which  was  approximately  uniform  over  the  cross-section, 
was  verified  using  PLIF.  No  liquid  droplets  were  observed  in  the  PLIF  plane  near  the  inlet  of  the 
test  section.  The  molecular  weight  of  the  jet  fuels  ranged  from  153.9  to  161  g/mol,  and  the 
hydrogen-to-carbon  ratio  of  the  jet  fuels  ranged  from  2. 017  to  2.176. 

Lean-blowout  measurements  were  conducted  on  four  types  of  stainless-steel  bluff-body 
flameholders:  v-gutters  with  diameters  of  0.95  -  3.8  cm,  a  cylinder,  a  plate,  and  a  square 
cylinder.  The  cylinder,  plate,  and  square  cylinder  were  3.8  cm  wide.  Blowout  was  achieved  by 
first  reaching  the  desired  inlet  conditions  and  then  decreasing  the  fuel  flow  until  the  flame  was 
very  near  blowout  but  stable.  The  fuel  flow  was  then  decreased  slowly  until  the  flame  blew  out, 
and  a  data  point  was  recorded. 

Correlation  Setup 

The  data  points  from  these  experiments  were  combined  with  data  points  obtained  from  the 
publications  of  other  authors,  including  Ballal  and  Lefebvre  [7],  DeZubay  [2],  King  [3],  Potter 
and  Wong  [18],  Williams  [1],  Zukoski  [4],  Chauduri  [14],  Yamaguchi  [19],  and  Barrere  [20], 
Table  1  provides  a  summary  of  the  experimental  parameters  from  each  study.  Hopane,  the  fuel 


30 

Approved  for  public  release;  distribution  unlimited. 


shown  for  the  DeZubay  work,  contained  95%  propane  and  5%  butane  and/or  ethane.  The  jet 
fuels  used  for  the  AFRL  dataset  were  Fuel  6169,  a  petroleum-derived  JP-8. 

Fuel  6308,  a  tallow-derived  HRJ  fuel,  and  Fuel  7721,  which  is  a  50-50  blend  of  6169  and  a 
camelina-derived  HRJ  fuel.  An  input  file  of  more  than  1100  individual  points  included  the 
relevant  parameters  from  each  study,  such  as  inlet  conditions  and  flameholder  characteristics. 

The  least-squares  curve-fit  function  in  the  Matlab  optimization  toolbox  was  used  to  find  an 
optimal  correlation.  The  desired  equation  form  was  selected  before  running  the  correlation  code. 

In  order  to  evaluate  how  various  parameters  increased  or  decreased  the  effectiveness  of  the 
correlation,  the  adjusted  R-squared  value  for  each  correlation  was  calculated  according  to  Eq.  15. 

R 2  SSj(n-p) 

SST/{n- 1)  (15) 

In  this  equation  SSe  is  the  sum  of  the  squares  of  the  residuals,  SSx  the  total  sum  of  the  squares  of 
the  original  dataset,  n  the  number  of  data  points,  and  p  the  number  of  parameters  used  in  the 
correlation. 

The  adjusted  R-squared  statistic  is  more  appropriate  than  the  traditional  R-squared  statistic  for 
comparing  correlations  that  contain  varying  numbers  of  parameters.  A  least-squares  curve-fit 
algorithm  will  run  until  the  residuals  have  been  reduced  as  much  as  possible.  If  a  parameter  is 
introduced  that  is  not  relevant  to  the  correlation,  the  curve  fit  will  simply  reduce  its  significance 
until  it  is  negligible,  and  the  change  to  the  goodness  of  fit  for  that  correlation  will  remain 
unchanged.  If  the  parameter  is  relevant,  the  optimization  process  will  assign  an  appropriate 
factor  such  that  the  goodness  of  fit  should  improve.  Because  the  least-squares  algorithm  will 
always  reduce  the  residuals  as  much  as  possible,  introducing  a  new  parameter  should  never 
worsen  the  fit,  and  will  therefore  never  decrease  the  traditional  R-squared  value.  The  adjusted 
R-squared  statistic  penalizes  added  parameters  to  correct  for  this  bias,  providing  a  better  statistic 
for  comparing  correlations  [21],  All  R-squared  statistics  reported  throughout  this  work  will  be 
the  adjusted  R-squared  values. 

In  addition  to  calculating  the  R-squared  statistic,  the  residual  for  each  point  was  plotted  against 
the  predicted  equivalence  ratio  at  lean  blowout  to  determine  whether  trends  were  present  within 
the  residuals.  If  trends  are  present,  it  indicates  that  the  model  is  inadequate,  possibly  due  to 
phenomena  that  are  not  fully  captured  in  the  correlation  [21],  For  the  correlations  in  this  study, 
the  residuals  were  randomly  distributed  over  the  range  of  predicted  equivalence  ratios. 

Cantera  [29]  was  used  to  calculate  the  laminar  flame  speeds  and  ignition  delay  times  for  selected 
data.  Two  hundred  data  points  over  a  range  of  pressures,  temperatures,  velocities,  Reynolds 
numbers,  and  flameholder  geometries  were  selected  from  the  blowout  experiments  on  propane 
and  methane  for  the  laminar  flame  speed  computation.  Propane  and  methane  were  chosen 
because  chemical  mechanisms  that  are  validated  at  relevant  conditions  for  these  fuels  are  readily 
available.  The  GRI-Mech  3.0  [22]  mechanism  was  used  for  the  methane  chemistry,  and  a 
mechanism  developed  by  Gokulakrishnan  et  al.  [23]  was  used  for  the  propane  chemistry.  GRI- 
Mech  3.0  is  an  optimized  mechanism  designed  to  model  natural-gas  combustion  with  53  species 
and  325  reactions  [22],  while  the  detailed  propane  mechanism  [23]  consists  of  136  species  and 
966  reactions.  The  propane  mechanism  has  been  validated  against  experimental  data  at 


31 

Approved  for  public  release;  distribution  unlimited. 


temperatures  pertinent  to  the  present  study  (300  -  650  K),  as  shown  in  Figure  4  [23].  Cantera 
was  employed  to  calculate  the  laminar  flame  speed  using  the  inlet  conditions  at  each  data  point. 
For  the  data  evaluated,  the  flame  speeds  ranged  from  6.73  to  55.1 1  cm/s.  These  data  were  then 
processed  through  the  Matlab  optimization  toolbox  in  a  manner  similar  to  that  described 
previously. 


Equivalence  Ratio 

Figure  16.  Validation  of  Cse  Propane  Mechanism  against  Flame  Speed  Experimental  Data  [23] 

Since  methane  and  propane  have  long  ignition  delay  times  at  atmospheric  pressure  and  relatively 
low  temperature  (below  600  K),  the  dataset  for  jet  fuels  was  selected  to  investigate  the  blowout 
correlation  with  ignition  delay  time.  A  detailed  surrogate  kinetic  mechanism  developed  by 
Gokulakrishnan  et  al.  [30]  was  used  to  compute  the  ignition  delay  time  of  jet  fuels.  This 
mechanism  uses  a  surrogate  mixture  for  each  jet  fuel  composed  of  n-dodecane,  n-decane,  iso¬ 
octane,  and  propyl-benzene.  The  ratio  of  the  various  components  is  determined  by  matching  the 
molecular  weight,  hydrogen-to-carbon  ratio,  derived  cetane  number,  and  threshold  sooting  index 
of  the  target  fuel.  Surrogates  that  mimic  these  global  fuel  properties  exhibit  chemical  kinetics- 
related  behavior  very  similar  to  that  of  the  target  fuels  they  are  meant  to  emulate  [24,  25],  The 
ignition  delay  time  was  calculated  for  each  point  based  on  its  inlet  conditions  using  Cantera  [29]. 
The  oxidizer  composition  was  calculated  using  a  Matlab  function  that  relies  on  combustion 
chemistry  and  thermodynamic  principles  of  equilibrium  to  calculate  the  chemical  make-up  of  the 
inlet  gases  after  vitiation.  The  Matlab  least-squares  curve-fit  function  then  used  the  resulting 
data  to  find  the  desired  correlation.  For  the  data  evaluated,  the  ignition  delay  times  ranged  from 
412  to  5728  ms. 

2.2.3  Correlation  Results 

In  a  previous  study  of  Huelskamp  et  al.  [16],  a  portion  of  the  current  dataset  was  used  to  explore 
which  flow  parameters  and  dimensionless  numbers  were  relevant  in  a  correlation  for  predicting 
lean  blowout.  The  smaller  dataset  did  not  account  for  vitiation  effects,  nor  did  it  include  the 
alternative  jet  fuel  data.  Because  this  work  found  the  most  relevant  parameters  to  be  lip  velocity, 
flameholder  diameter,  pressure,  and  temperature,  these  were  the  initial  factors  considered  in  the 
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present  analysis.  In  addition,  it  was  found  that  employing  a  ratio  of  lip  velocity  to  flameholder 
diameter  instead  of  using  each  factor  individually  was  only  slightly  detrimental  to  the  goodness 
of  fit  of  the  correlation,  and  this  ratio  is  recognized  as  a  fluid-mechanic  timescale.  As  such,  a 
ratio  of  the  two  values  is  used  in  the  following  correlations.  In  these  correlations,  velocity  (U)  is 
in  m/s,  pressure  (P)  in  atm,  diameter  (D)  in  m,  and  temperature  (T)  in  K. 

Figure  17  shows  the  correlation  results  with  all  of  the  available  data.  The  data  are  plotted  with 
the  predicted  equivalence  ratio  (or  phi)  at  lean  blowout  (LBO)  on  the  x-axis  and  the  actual 
equivalence  ratio  at  lean  blowout  on  the  y-axis.  The  black  line  running  through  the  plot 
represents  a  “perfect”  correlation  (i.e.,  the  predicted  phi  at  LBO  is  equal  to  the  experimental  phi 
at  LBO).  The  plot  legend  indicates  the  fuel  used  in  each  experiment.  Equation  16  is  the 
correlation  equation.  The  R-squared  value  from  this  correlation  is  0.852. 
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Figure  16.  Correlation  Results  with  D/U,  Temperature,  and  Pressure  as  Factors  (R2=0.852) 
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Despite  the  relatively  good  fit  of  the  data  to  the  correlation  equation,  some  effects  are  visibly  not 
captured.  Close  examination  of  Figure  5  reveals  that  data  points  obtained  from  experiments 
using  fuels  other  than  propane  appear  either  to  fall  away  from  the  “perfect- fit”  black  line,  as  in 
the  case  of  methane  and  alternative  jet  fuel  data,  or  to  follow  a  slightly  different  slope  on  the 
plot,  as  in  the  case  of  the  JP-4  data.  This  observation  led  to  exploring  the  addition  of  fuel 
characteristics  to  the  correlation. 


As  discussed  in  the  Experimental  Setup  Section,  researchers  have  successfully  developed 
surrogate  jet  fuel  formulations  based  on  molecular  weight,  hydrogen-to-carbon  ratio,  threshold 
sooting  index  (TSI),  and  derived  cetane  number  (DCN)  [24,  25].  Of  these  parameters,  both 
molecular  weight  and  hydrogen-to-carbon  ratio  are  easily  obtained  for  all  of  the  fuels  in  the 
dataset  for  this  study;  therefore,  the  molecular  weight  (MW)  and  hydrogen-to-carbon  ratios 
(H/C)  for  each  fuel  were  added  to  the  correlation.  The  parameters  were  added  separately  so  that 
the  effect  of  each  could  be  discerned. 

Figure  18  displays  the  data  plotted  after  adding  the  fuel  molecular  weight  to  the  correlation 
equation,  as  shown  in  Eq  17. 
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Figure  17.  Correlation  Results  with  Molecular  Weight  as  Fuel  Factor  (R2=0.858) 
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The  R-squared  value  from  this  correlation  is  slightly  higher  than  that  from  the  correlation 
without  molecular  weight,  increasing  the  value  by  0.006.  The  exponent  on  molecular  weight  in 
the  correlation  is  also  quite  small.  Examination  of  Figure  6  shows  that  the  only  readily  apparent 
benefit  of  the  additional  parameter  is  to  bring  the  Zukoski  methane  data  [4]  closer  to  the  perfect 
fit-line.  The  Zukoski  gasoline  data  [4]  and  the  AFRL  jet  fuel  data  actually  appear  to  drift  farther 
from  the  perfect- fit  line  with  the  addition  of  molecular  weight  in  the  correlation. 

Figure  19  shows  the  data  plotted  with  the  addition  of  the  hydrogen-to-carbon  ratio  of  each  fuel  to 
the  correlation  equation,  shown  in  Eq.  18.  The  R-squared  value  of  the  correlation  is  0.861, 
which  is  higher  than  the  R-squared  value  of  the  correlation  with  molecular  weight  and  the  R- 
squared  value  of  the  initial  correlation,  0.852.  The  addition  of  the  hydrogen-to-carbon  ratio 
appears  to  bring  the  methane  data  closer  to  the  perfect-fit  line  and  does  not  have  the  negative 
influence  that  the  molecular  weight  had  on  the  AFRL  jet  fuel  data. 
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Figure  18.  Correlation  Results  with  Hydrogen-to-Carbon  Ratio  as  Fuel  Factor  (R2=0.861) 
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Although  the  addition  of  the  hydrogen-to-carbon  ratio  improves  the  R-squared  value  and  brings 
the  methane  data  closer  to  the  perfect- fit  line,  the  AFRL  jet  fuel  data  continue  to  diverge  from 
the  perfect  fit.  The  AFRL  jet  fuel  data  as  well  as  King’s  JP-4  data  [3]  were  obtained  in  vitiated 
flows.  Vitiation  introduces  a  variety  of  chemical  species  that  can  affect  the  combustion  reaction 
and,  in  turn,  the  Damkdhler  number  and  the  lean-blowout  process.  In  addition,  vitiation  reduces 
the  amount  of  oxygen  available  in  the  combustion  reaction.  Researchers  have  found  that 
vitiation  results  in  higher  equivalence  ratios  at  lean  blowout  because  of  decreased  reactivity  from 
reduced  oxygen  levels  [26]. 

In  an  attempt  to  capture  these  effects,  the  vitiation  level  was  included  in  the  correlation  by 
calculating  the  percentage  of  oxygen  in  the  oxidizing  gas  using  the  Matlab  combustion  chemistry 
code  discussed  previously.  In  the  case  of  the  AFRL  data,  this  calculation  could  be  done  directly 
using  actual  air  and  fuel  flows  recorded  by  the  combustion  facility.  The  oxygen  level  for  King’s 
data  was  obtained  by  varying  the  vitiator  fuel-air  ratio  in  the  code  until  the  gas  temperature 
matched  that  of  the  inlet  temperature  reported  by  King  [3].  The  correlation  was  then  run  again 
with  both  the  oxygen  level  and  the  hydrogen-to-carbon  ratio  included.  Figure  20  shows  the  data 
plotted  with  the  new  correlation;  the  correlation  equation  is  shown  in  Eq.  19.  The  R-squared 
value  of  the  correlation  has  improved  to  0.873,  from  a  previous  value  of  0.861.  The  AFRL  jet 
fuel  data  has  moved  closer  to  the  perfect-fit  line. 
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Figure  20.  Correlation  Results  with  H/C  and  Oxygen  Content,  (R2=0.873) 
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2.2.3. 1  Correlation  with  Laminar  Flame  Speed 

Cantera  was  used  to  calculate  the  laminar  flame  speeds  for  200  data  points  taken  on  propane  and 
methane  over  a  range  of  pressures,  temperatures,  velocities,  Reynolds  numbers,  and  flameholder 
geometries.  These  data  were  then  processed  with  the  Matlab  optimization  toolbox  to  investigate 
the  effect  of  laminar  flame  speed  on  blowout.  Adjusted  R-squared  values  were  also  calculated. 

As  discussed  previously,  Kariuki  et  al.  [10]  proposed  the  use  of  the  parameter  in  Eq.  20  to 
predict  when  lean  blowout  will  occur.  This  parameter  is  derived  from  a  rearrangement  of  the 
turbulent  premixed  flame  extinction  theory  originally  proposed  by  Radhakrishnan  et  al.  [9],  In 
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Eq  (8),  Ub  is  the  bulk  velocity,  d  the  flameholder  diameter,  Sl  the  laminar  flame  speed,  and  v  the 
kinematic  viscosity  of  the  reactants. 
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Kariuki  et  al.  [10]  plotted  this  value  against  the  flame  power  for  swirl-stabilized  flames  and 
found  that  the  values  for  premixed  flames  were  near  unity.  The  dataset  for  which  laminar  flame 
speed  was  calculated  in  the  current  study  was  used  in  an  attempt  to  verify  the  results  of  previous 
researchers.  Figure  21  shows  the  parameter  of  Kariuki  et  al.  plotted  against  the  equivalence  ratio 
at  lean  blowout.  While  the  propane  data  do  indeed  produce  a  parameter  near  unity — especially 
at  an  equivalence  ratio  over  0.7 — the  methane  data  resulted  in  a  parameter  value  of  up  to  10. 

The  data  were  then  run  through  the  Matlab  correlation  code.  To  establish  a  baseline,  the  initial 
parameters  used  were  D/U,  pressure,  and  temperature.  Equation  2 1  shows  the  resulting 
correlation  equation,  and  the  data  are  plotted  in  Figure  22.  The  R-squared  value  was  0.802.  The 
methane  data  in  Figure  10  clearly  do  not  follow  the  same  trend  as  the  propane  data.  In  addition, 
it  appears  that  there  are  two  distinct  curves  within  the  methane  data,  each  corresponding  to  a 
different  D/U  ratio. 
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Figure  19.  Kariuki  Parameter  Plotted  against  Equivalence  Ratio  at  Lean  Blowout 


Figure  20.  Flame  Speed  Dataset  Correlation  Results  with  D/U,  Temperature,  and  Pressure  as 

Factors  (R2=0.802) 
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Because  laminar  flame  speed  has  been  indicated  by  previous  researchers  [9,  10]  as  a  critical 
parameter  in  blowout,  it  was  used  to  replace  the  pressure  and  temperature  in  the  correlation  as 
the  chemical  “timescale.”  The  resulting  correlation  equation  is  Eq  22.  Including  laminar  flame 
speed  improved  the  R-squared  value  from  0.802  to  0.846. 
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Equation  22  reveals  multiple  unexpected  outcomes.  The  fluid-dynamic  timescale  appears  to 
have  lost  its  significance  in  the  blowout  process.  In  addition,  the  laminar  flame  speed  is  in  the 
numerator,  suggesting  that  the  higher  the  laminar  flame  speed,  the  less  stable  the  flame — which 
is  in  opposition  to  most  physical  explanations  of  flame  stability.  A  more  likely  explanation 
emphasizes  the  need  for  caution  when  using  correlations  to  explore  physical  phenomena  in  large 
datasets. 

It  is  very  probable  that  the  correlation  is  not  actually  predicting  lean  blowout  but  rather  is 
creating  a  laminar  flame  speed  curve.  Laminar  flame  speed  increases  with  equivalence  ratio 
when  the  flame  is  lean,  as  demonstrated  previously  in  Figure  17.  The  correlation  algorithm 
recognized  this  and  simply  correlated  the  flame  speed  to  the  equivalence  ratio.  If  more  data  were 
available  where  the  flame  speeds  differed  substantially  at  the  same  equivalence  ratio,  the 
equivalence  ratio/flame  speed  link  might  be  broken,  and  the  correlation  could  be  run  to  predict 
lean  blowout.  With  the  current  dataset,  the  laminar  flame  speed  cannot  be  used  to  predict 
blowout. 

2.2.Z.2  Correlation  with  Ignition  Delay  Time 

Ignition  delay  time  is  another  chemical  timescale  that  has  been  proposed  as  being  critical  to  the 
blowout  process;  therefore,  a  subset  of  the  data  was  used  to  correlate  with  ignition  delay  time. 
The  methane  and  propane  dataset  used  for  the  laminar  flame  speed  calculations  could  not  be 
used  to  study  the  effect  of  ignition  delay  time  because  of  the  extremely  long  ignition  delay  times 
at  atmospheric  pressure  and  relatively  low  temperatures  (below  600  K).  Instead,  the  dataset 
chosen  for  the  ignition  delay  time  correlation  was  composed  of  data  taken  on  a  petroleum- 
derived  JP-8  fuel  and  a  tallow-derived  HRJ  fuel,  which  have  relatively  shorter  ignition  delay 
times.  The  chemical  mechanism  employed  [30]  uses  surrogate  mixtures  to  model  each  fuel 
using  varying  ratios  of  n-dodecane,  n-decane,  iso-octane,  and  propyl-benzene.  The  ignition 
delay  time  for  each  of  the  40  points  in  the  dataset  was  calculated  based  on  its  inlet  conditions 
using  Cantera  [29],  Because  each  of  these  points  was  taken  on  vitiated  air,  the  oxidizer 
composition  was  calculated  using  a  Matlab  function  that  relies  on  combustion  chemistry  and 
thermodynamic  principles  of  equilibrium  to  calculate  the  chemical  make-up  of  the  inlet  gases 
after  vitiation.  The  Matlab  least-squares  curve-fit  function  was  then  used  to  find  the  desired 
correlation. 

The  correlation  was  first  run  using  only  the  inlet  conditions  to  establish  a  baseline  comparison. 
The  data  correlate  quite  well  with  only  the  inlet  conditions,  probably  because  less  variation  was 
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present  in  the  dataset  than  would  be  the  case  with  more  data  points.  The  pressure  varied  by  only 
around  0.04  atm,  and  the  temperature  ranged  from  750  to  910  K.  In  addition,  all  points  were 
obtained  at  the  AFRL  High  Pressure  Combustion  Research  Facility,  which  eliminated  some 
variability  due  to  facility-related  inconsistencies.  The  correlation  data  are  plotted  in  Figure  23, 
and  the  correlation  equation  is  Eq  23.  The  R-squared  value  is  0.942. 
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Figure  21.  Ignition  Delay  Dataset  Correlation  Results  with  D/U,  Temperature,  and  Pressure  as 

Factors  (R2=0.942) 

Because  both  pressure  and  temperature  affect  ignition  delay  time,  they  were  replaced  in  the  next 
correlation  with  the  ignition  delay  time  (IDT)  as  shown  in  Eq  24.  The  data  are  plotted  in  Figure 
24.  When  the  ignition  delay  time  was  substituted  into  the  correlation,  the  R-squared  value 
decreased  to  0.918.  However,  the  data  continues  to  fall  very  near  the  perfect- fit  line. 


Figure  22.  Ignition  Delay  Dataset  Correlation  Results  with  D/U  and  IDT  as  Factors  (R2=0.918) 
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To  verify  that  the  data  were  not  solely  dependent  on  the  fluid-mechanic  timescale,  the  correlation 
was  also  run  with  the  D/U  parameter  only,  as  shown  in  Eq  25.  This  resulted  in  an  R-squared 
value  of  0.381.  In  addition,  the  data  no  longer  fell  as  close  to  the  perfect-fit  line,  as  shown  in 
Figure  25. 
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Figure  23.  Ignition  Delay  Dataset  Correlation  Results  with  D/U  as  the  Only  Factor  (R2=0.381) 
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The  data  are  more  scattered  when  no  chemical  timescale  parameters  are  included.  Although  the 
pressure  and  temperature  are  adequate  to  fit  the  data,  they  are  not  direct  timescales  that  can  be 
included  in  a  Damkohler-type  number.  While  ignition  delay  time  appears  to  be  an  adequate 
substitute  for  this  particular  dataset,  the  dataset  was  relatively  small.  All  points  were  taken  in  the 
High  Pressure  Combustion  Research  Facility  at  AFRF  within  a  fairly  small  range  of  inlet 
conditions  and  fuel  properties.  Because  the  data  lack  a  lot  of  variation,  the  correlation  does  a 
very  good  job  of  predicting  blowout  without  specifically  identifying  a  chemical  timescale; 
therefore,  it  is  difficult  to  discern  whether  including  ignition  delay  time  would  have  a  more  or 
less  significant  impact  in  a  larger  dataset.  However,  for  the  dataset  used  in  this  study,  the 
ignition  delay  time  did  provide  an  adequate  representation  of  chemical  timescale. 

2.2.3.3  Correlation  Validation 

To  assess  the  accuracy  of  the  correlation,  the  correlation  including  hydro  gen- to-carbon  ratio  and 
oxygen  content  was  used  to  predict  the  blowout  of  bluff-body  flames  at  an  additional  47 
conditions.  The  blowout  points  were  selected  from  propane  and  jet  fuel  experiments  run  in  the 
Atmospheric  Pressure  Combustion  Faboratory  at  AFRF  as  well  as  from  propane  experiments 
conducted  by  Yamaguchi  et  al.  [19].  Figure  14  displays  a  plot  of  the  validation  points,  shown  in 
yellow,  as  well  as  the  points  used  in  the  correlation,  the  equation  of  which  is  shown  again  in  Eq 
26. 
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The  average  error  of  the  prediction  within  the  validation  dataset  was  10.57%.  The  AFRL 
propane  and  jet  fuel  datasets  lie  near  the  perfect-fit  line  and  have  an  average  error  of  6.92%  and 
5.08%,  respectively. 


The  data  of  Yamaguchi  et  al.  also  lie  relatively  near  the  perfect-fit  line,  with  the  exception  of  an 
obvious  outlier  at  an  actual  blowout  of  0.4.  This  point  was  also  something  of  an  outlier  in  the 
original  data  of  Yamaguchi  et  al  [19].  Of  particular  interest  concerning  the  Yamaguchi  data  is 
that  flameholder  configurations  employing  multiple  bluff  bodies  were  used.  The  blockage  ratio 
of  the  experiment  was  kept  constant  by  changing  the  flameholder  diameter  when  the  number  of 
flameholders  was  changed  [19].  When  calculating  the  predicted  equivalence  ratio  at  lean 
blowout,  only  the  actual  diameter  of  the  individual  flameholders  was  taken  into  account.  The 
correlation  predicted  the  lean  blowout  of  the  Yamaguchi  points  with  an  average  error  of  19.54%. 
If  the  outlier  is  excluded,  the  average  error  falls  to  16.78%. 
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Figure  24.  Correlation  Results  with  H/C  and  Oxygen  Content,  Showing  Validation  Points 


The  average  error  in  the  predicted  equivalence  ratio  for  all  points,  including  those  used  to  form 
the  correlation  equations,  was  8.10%.  The  maximum  error  was  64.00%,  which  was  the  error  on 
a  point  taken  on  methane.  Methane  points  were  generally  associated  with  higher  errors,  with  an 
average  of  19.88%. 

2.2.4  Discussion 


Several  of  the  parameters  in  the  correlation  shown  in  Eq  26  are  logical  extensions  of  the 
Damkohler  number  theories  presented  in  previously  published  literature.  D/U  is  the  fluid 
mechanic  timescale,  where  U  is  the  lip  velocity  and  D  the  flameholder  diameter.  This  finding 
supports  the  conclusions  of  many  past  researchers  [5,  8,  12]  and  implies  that  the  mixing  time  in 
the  shear  layer  between  the  hot  products  in  the  recirculation  zone  and  the  fresh  reactants  flowing 
past  the  flameholder  is  an  important  timescale.  As  the  fresh  reactants  flow  past  the  recirculation 
zone,  they  mix  in  the  shear  layer  with  hot  combustion  products.  If  the  reactants  do  not  mix  and 
react  sufficiently  in  the  amount  of  time  that  they  spend  in  the  shear  layer,  the  flame  cannot 
propagate.  Because  the  dimensions  of  the  recirculation  zone  generally  scale  with  the  width  of 
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the  flameholder,  the  width  of  the  flameholder  can  be  substituted  for  the  recirculation  zone  length, 
although  it  introduces  some  variability  to  the  correlation. 

Both  pressure  and  temperature  contribute  to  the  chemical  timescale  in  the  Damkohler  number. 

In  the  correlation  equation  shown  in  Eq.  26,  the  equivalence  ratio  at  lean  blowout  decreases  as 
the  pressure  and  temperature  increase.  In  general,  as  the  pressure  and  temperature  of  a 
combustible  mixture  increase,  the  reactivity  of  that  mixture  increases.  If  the  reactivity  of  the 
mixture  in  a  bluff-body  flame  increases,  the  chemical  timescale  should  decrease,  resulting  in  a 
flame  that  should  be  able  to  propagate  in  more  adverse  flow  conditions. 

For  this  dataset,  the  molecular  weight  appears  to  have  a  minimal  effect  on  the  equivalence  ratio 
at  lean  blowout.  However,  the  hydrogen-to-carbon  ratio  does  improve  the  correlation.  The 
hydrogen-to-carbon  ratio  is,  to  some  extent,  an  indication  of  the  tendency  of  the  fuel  to  ignite,  or 
its  cetane  number.  The  cetane  number  of  chemical  components  typically  present  in  jet  fuels 
shows  the  following  trend:  normal  alkanes  >  branched  alkanes  >  cycloalkanes  >  aromatics  [27], 
As  demonstrated  in  Table  4,  the  hydrogen-to-carbon  ratio  of  these  structures  follows  a  similar 
trend  and,  therefore,  is  an  indicator  of  the  relative  ignition  time  of  the  fuel.  If  the  ignition  delay 
time  is  an  adequate  representation  of  the  critical  chemical  timescale  in  the  Damkohler  number, 
then  it  logically  follows  that  the  cetane  number — and,  therefore,  the  hydrogen-to-carbon  ratio — 
would  have  an  influence  in  the  correlation  for  jet  fuels. 

Table  4.  Hydrocarbon  Structures  and  Their  Relative  Cetane  Numbers  and  H/C  Ratios 


<D 

n 

E 

D 

z 

Q) 

C 

re 

*-* 

a> 

O 

U) 

£ 

'</> 

re 

a> 

»_ 

o 

c 


A  relationship  between  cetane  number  and  lean  blowout  was  also  found  by  Colket  et  al.  [28] 
who  used  experimental  combustor  data  taken  with  a  variety  of  jet  fuels.  They  found  that  fuels 
with  higher  cetane  numbers  blew  out  at  lower  equivalence  ratios,  as  shown  in  Figure  27.  In 
addition,  they  found  that  fuels  with  more  n-alkanes,  which  have  high  hydrogen-to-carbon  ratios, 
are  less  likely  to  blow  out. 
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Figure  25.  Derived  Cetane  Number  against  Equivalence  Ratio  at  LBO  [28] 


In  the  correlation  shown  in  Eq  14,  the  hydro  gen- to-carbon  ratio  is  in  the  denominator,  which 
implies  that  the  higher  the  hydrogen-to-carbon  ratio,  the  more  resistant  the  flame  is  to  lean 
blowout.  In  jet  fuels,  a  higher  hydrogen-to-carbon  ratio  indicates  more  n-alkanes  and  a  higher 
cetane  number;  therefore,  a  more  stable  flame  would  indeed  be  anticipated. 

The  last  parameter  included  in  Eq  26  is  the  oxygen  concentration  of  the  incoming  air  after 
vitiation.  Less  oxygen  content  in  the  air  should  lengthen  the  chemical  timescales  involved  and, 
thereby,  increase  the  equivalence  ratio  at  lean  blowout.  This  behavior  is  confirmed  by  the 
presence  of  the  oxygen  content  in  the  denominator  of  Eq  26. 


2.2.5  Conclusions 


Lean  blowout  data  was  collected  on  propane  and  jet  fuel  bluff-body  stabilized  flames  and  was 
combined  with  data  taken  from  past  literature  to  create  a  database  of  over  1 100  data  points. 
Cantera  was  used  to  calculate  the  laminar  flame  speeds  for  200  data  points  taken  on  propane  and 
methane  over  a  range  of  pressures,  temperatures,  velocities,  Reynolds  numbers,  and  flameholder 
geometries.  Ignition  delay  times  were  also  calculated  for  40  data  points  taken  on  conventional 
and  alternative  jet  fuels.  A  least-squares  curve-fit  algorithm  was  run  using  Matlab  to  find  a 
correlation  that  could  accurately  predict  the  equivalence  ratio  at  lean  blowout  for  each  dataset. 

The  best  correlation  for  the  overall  dataset  included  pressure,  temperature,  the  ratio  of 
flameholder  diameter  to  lip  velocity,  oxygen  level,  and  the  hydrogen-to-carbon  ratio  of  the  fuel 
as  factors.  The  R-squared  value  of  this  correlation  was  0.873.  The  exponents  on  the  factors 
indicated  that  high  pressure  and  temperature  lowered  the  equivalence  ratio  at  lean  blowout,  as 
did  high  levels  of  oxygen.  These  factors  contribute  to  the  chemical  timescale  in  the  Damkohler 
number.  The  correlation  also  indicated  that  the  equivalence  ratio  at  lean  blowout  decreases  as 
the  hydrogen-to-carbon  ratio  increases.  In  jet  fuels,  high  hydrogen-to-carbon  ratios  indicate 
increased  levels  of  alkanes,  which  have  higher  cetane  numbers  and  shorter  ignition  delay  times. 
The  ratio  of  flameholder  diameter  to  lip  velocity  is  recognized  as  the  fluid  mechanic  timescale  in 
the  Damkohler  number. 

The  correlation  for  the  ignition  delay  dataset  with  pressure,  temperature,  and  the  ratio  of 
flameholder  diameter  to  lip  velocity  as  factors  had  an  R-squared  value  of  0.942.  This  relatively 
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small,  uniform  dataset  correlates  very  well.  When  ignition  delay  time  replaces  pressure  and 
temperature,  the  data  continues  to  correlate  well,  with  an  R-squared  value  of  0.918.  For  this 
limited  dataset,  the  ignition  delay  time  is  an  adequate  representation  of  the  chemical  timescale. 
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